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Visual  observations  made  on  a  pulsed  toroidal  plasma  indicate  that 
the  current  carrying  plasma  sheet  breaks  up  into  filaments  parallel  to  the 
direction  of  the  current  and  perpendicular  to  the  direction  of  the  magnetic 
field.  This  process  involves  the  annihilation  and  fusion  of  initially  opposite 
magnetic  flux  lines  and  leads  to  a  final  configuration  which  corresponds  to 
a  lower  energy  state.  Such  a  state,  however,  would  not  be  accessible  in  the 
absence  of  dissipation  since  its  attainment  requires  slippage  of  the  plasma 
across  magnetic  flux  lines.  The  operating  mechanism  must  then  be 
associated  with  a  resistive  instability  of  the  type  first  predicted  by  Furth, 
Killeen,  and  Rosenbluth  and  also  by  Coppi,  Greene  and  Johnson  in  Physics 
of  Fluids,  Volume  6,  1963.  Observations  of  resistive  instabilities  under 
laboratory  conditions  is  of  particular  interest  because  these  modes  are  held 
responsible  for  many  astrophysical  and  geophysical  phenomena  such  as  solar 
flares  and  auroras  as  well  an  anomalous  diffusion  rates  in  fusion  devices. 

The  experimental  setup  consists  of  a  plasma  filled  ring-shaped 
container  of  rectangular  cross  section,  where  a  discharge  is  induced  by 
means  of  an  iron  cored  transformer.  The  time  development  of  the  current 
distribution  in  the  plasma  is  determined  with  the  help  of  a  combination  of 
magnetic  probes  and  an  S.T.L.  Image  Converter  Camera.  Optical  spectro¬ 
scopy  is  used  to  calculate  the  temporal  variation  of  electron  number  density 
and  temperature. 


iii 


f 


It  is  found  that  the  plasma  forms  near  the  inner  cylindrical  wall, 
where  the  primary  winding  is  located,  and  moves  radially  outward  as  time 
progresses.  The  radial  speed  depends  on  the  mass  density  and  magnetic 
field  associated  with  the  plasma  in  accordance  with  the  snowplow  model. 
Concurrently  with  the  outward  motion  a  resistive  instability  develops,  which 
breaks  the  current  sheet  into  separate  rings.  The  wavelength  of  the 
instability  and  therefore  the  number  of  rings  depend  on  the  geometry  of  the 
primary  winding  and  on  the  initial  perturbation  it  introduces  in  the  field 
distribution.  It  was  found  that  for  a  primary  winding  having  N  equal  segments, 
the  initially  uniform  plasma  sheet  split  into  N  +  1  rings. 

Theory  predicts  that  the  "tearing"  and  "rippling"  modes  of  a 
resistive  instability  develop  with  an  exponential  growth  rate  that  is  a 
_ function  of  the  resistivity,  mass  density,  magnetic  flux  density,  and  plasma 
thickness.  The  experiments  confirm  this  functional  dependence.  Also 
consistent  with  the  theory  is  the  influence  of  wavenumber  on  the  growth  rate 
for  the  two  modes  mentioned  above. 
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CHAPTER  1 


INTRODUCTION 

There  are  many  reasons  for  investigating  instabilities  and  the 
confinement  of  plasmas.  These  problems  are  basic  to  various  physical 
phenomena,  proposed  fusion  devices,  propulsion  systems  and  MHP 
gen»iators.  At  present,  there  is  much  interest  in  resistive  instabilities 
since  some  plasma  configurations  which  may  be  stable  in  the  absence  of 
dissipation  become  instable  for  finite  electrical  resistivities.  In  particular, 
instabilities  of  the  resistive  type  are  held  responsible  for  astrophysical 
and  geophysical  phenomena  such  as  solar  flares  and  auroras  as  well  as 
anomalous  diffusion  rates  in  fusion  devices. 

The  first  theoretical  papers  written  on  the  topic  of  resistive 

instabilities  where  written  by  Furth,  Killeen,  and  Rocenbluth*,  and  Coppi, 
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Greene,  and  Johnson  in  Experimental  evidence  of  a  resistive 

instability  in  a  theta  pinch  was  reported  by  Bodin^  about  the  same  time. 

Bodin's  experiment  showed  the  plasma  breaking  up  into  several  filaments. 

4  5,  h 

Experimental  and  theoretical  papers  by  Kvartskava  and  Koma-ov 
respectively  are  interesting  because  they  investigate  spatially  periodic 
plasma  structures  which  can  arise  from  resistive  instabilities. 

It  is  recognized  that  there  should  always  be  a  balance  of  theoretical 
predictions  and  supporting  expei  imental  data  if  possible.  It  is  often  difficult, 
to  design  an  experiment  to  isolate  and  examine  a  specific  instability,  how¬ 
ever.  So  at  the  present  time,  material  devoted  to  the  theory  of  resistive 

instabilities  fa-  outweighs  experimental  evidence.  As  a  matter  of  fact  a 
7 

policy  paper  by  the  AEC  pointed  this  out  and  suggested  that  experimental 
efforts  be  intensified  towards  the  isolation  of  types  of  instabilities.  It  is 
with  this  viewpoint  the  thesis  was  conducted.  This  thesis  is  essentially 
experimental  and  presents  mainly  data  supporting  the  existence  of  resistive 
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instabilities  in  a  torotaal  piasma  composed  of  hydrogen  at  a  pressure 
of  500  microns.  Ai  a  much  lower  pressure,  50  microns,  data  supporting 
the  e  firt'-nce  of  Rayle igh- Taylor  insta oilities  is  presented. 

Observations  made  on  a  pulsed  toroidal  discharge  at  500  microns 
pressure  indicate  that  the  current  carrying  plasma  sheet,  see  Figure  I- 1 , 
breaks  up  into  filaments  parallel  tc  the  direction  of  the  current  and  perpen¬ 
dicular  to  the  direction  of  the  magnetic  field.  This  process  involves  the 
annihilation  and  fusion  of  initially  opposite  magnetic  flux  lines  and  lends  to 
a  final  configuration  which  corresponds  to  a  lower  energy  state.  Such  a  state, 
however,  would  not  be  accessible  in  the  absei  ~s  of  dissipation  since  its 
attainment  requires  slippage  of  the  plasma  across  magnetic  flux  lines.  Thus 
the  operating  mechanism  must  be  associated  with  a  resistive  instability. 

At  a  much  lower  pressure,  i.  e.  50  microns,  the  plasma  moves 
outward  at  a  high  radial  speed  in  accordance  with  the  snow  plow  model. 

The  acceleration  of  the  interface  between  the  magnetic  field  and  the  plasma 
produces  flutes  along;  the  surface  of  the  plasma.  See  Figure  1-2.  These 
fiutes  may  be  interpreted  as  Rayleigh- Taylor  instabilities,  and  the  wave- 

g 

length  of  the  flutes  can  be  related  to  viscous  damping 

Chapter  II  presents  some  physical  models  that  provide  insight  into  the 
theory  of  resistive  instabilities.  A  derivation  of  Rayleigh- Taylor  instabilities 
for  zero  resistivity  and  ot  the  damping  effect  produced  by  finite  viscosity 
follows.  A  description  of  the  torridal  discharge  and  method  of  excitation  is 
given  in  Chapter  III.  Chapter  IV  is  concerned  with  instrumentation  of  the 
experiment  which  includes  voltage  and  current  measurements,  magnetic 
probes,  high  speed  photography,  spectrograph,  and  scann ing  monochromate r. 
The  experimental  findings  are  given  in  Chapter  V.  These  include  data  from 
measurements  c'  magnetic  field,  electron  temperature,  resistivity,  electron 
density,  and  photographs  from  wHch  the  growth  rate  o,  the  various 
instabilities  may  bt  calculated. 


3 

Chapter  VI  presents  conclusions  about  the  measurements  of  various 
plasma  parameters,  the  observations  of  both  resistive  and  Rayleigh- Taylor 
instabilities  and  their  relationship  to  respective  theoretical  predictions. 


4 


MAGNETIC  FIELD 


PLASMA  SHEET  BEFORE  BREAK-UP 


MAGNETIC  FIELD 


PLASMA  SHEET  AFTER  BREAK-UP 
The  effect  of  a  resistive  instability 
F  igure  I- 1 


PLASMA 


■MAGNETIC  FIELD 
•INNER  WALL  OF  TOROID 

PLASMA  SHEET  MOVING  RACIALLY  OUTWARD 


DEVELOPMENT  OF  FLUTES 


The  effect  of  a  Rayleigh-Taylor 
Instability 

P'igure  1-2 


t 


6 


CHAPTER  II 


THEORETICAL  BEHAVIOUR  OF  THE  INSTABILITIES 


The  experimental  configuration  is  a  toroidal  current  carrying  plasma 
whose  cross-section  is  nearly  rectangular  or  sheet-like.  This  configuration 
is  described  in  Chapter  HI.  In  any  case  the  experimental  plasma  when  approx 
imated  as  a  t,heet  pinch  takes  on  the  geometry  as  discussed  by  Furth,  Kille  n 
and  Rosenbluth^  .  The  following  physical  interpretation  of  resistive  instabil¬ 
ities  follows  their  treatment. 

Shown  first  in  Figure  II- 1  is  a  sketch  of  a  sheet  pinch  indicating  the 
principal  axis  and  magnetic  field  supported  by  the  current  carrying  plasma. 
The  plasma  sheet  has  a  thickness  characterized  by  a. 

We  start  with  Ohm's  law,  which  is  : 

1J  =  E  +  VxB  2.1) 

and  if  we  choose  a  frame  of  reference  in  which  the  electric  field  intensity 
E  vanishes  (Bis  time  invarient)  equation  2.  1)  becomes 

T)  J  =  V  x  _B  2.2) 

The  motor  force,  xB,on-n  element  of  f  luid  is 

Fm  =  J  x  B  =  B  (  V-B)  -  V  R2  2.  3) 

It  is  seen  from  Figure  II- 1  and  equation  2.  3)  that  only  fluid  motion 
in  ihe  y  direction  results  in  a.  force  and  fur  ,he  rrr.ore  this  force  is  a  restor¬ 
ing  one.  So  equation  2.  3)  becomes 

VB^  ^n  Vv(B')2(ra)2 
Tl  ~  r\ 


Fm 


2.4) 


r 
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and  ca  is  tile  thickness  of  a  small  layer  about  B=0  where  B~B,(e:a), 
and  the  prime  -  ~  .Inspection  of  the  equation  above  shows  that  the  plasma 
motion  along  the  y  axis,  transverse  to  B,  is  hampered  .  In  fact  for  , 
the  plasma  is  frozen"  to  the  flux  lines.  On  the  other  hand,  at  a  point  where 
B  is  almost  zero  we  find  that  Fm  is  correspondingly  weak  so  that  the  plasma 
can  diffuse  easily.  The  region  in  which  this  happens  we  characterize  by  a 
length  of  ea  units. 


For  an  instability  to  occur  then,  a  force  must  be  exerted  of  magnitude 
|Fm|  and  opposite  in  direction  to  the  restoring  force.  Also  we  expect  that 
the  instability  starts  in  the  region  ea  .  The  power  associated  with  the  restor¬ 
ing  force  is  of  the  order 


P~  -  V- Fm 

2,2  2 

V  (BT  (ea)4 

or  P  ~  -1 -  2.5) 


Any  instability  that  occurs  will  be  observed  in  the  xy  plane.  In 
general,  the  wavelength  of  the  instability  (X  =  “  )  is  much  larger  than 
ea  so  that  the  direction  of  propagation  of  the  instability  is  along  the  x  axis 
rather  than  the  y  axis.  Therefore  the  instability  is  assumed  to  have  the 
following  functional  dependence  . 

F(x,y,t)  =  Fj (y)  e^  eXt  2.6) 

The  power  associated  with  the  instability  is  then 


P 


instab. 


_d_ 

dt 


2.  7) 


where 


V 


_JL_ 


ea 


since  a.  low  velocities  the  plasma  is  practically  incompressible. 
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Simplifying  equation  2.  7)  and  equating  it  to  the  restoring  power  v.e  have 


.  2  ,,2  „2  /«'»2/ 

4n  P  Vy  _  (B  )  (ea) 


k2(ea)2 


Then  the  thickness  of  the  unstable  layer  ea  is  given  by  the  following  : 


ea  ~ 


4n  nipn 

L  k2(B')2 


y« 


2.  8) 


In  this  experiment  two  specific  modes  connected  with  resistive 
inatabilities  are  of  interest.  The  modes  have  been  termed  the  "rippling" 
mode  and  the  "tearing"  mode  by  F.K.  R.  The  rippling  mode  has  been 
associated  with  the  gradient  in  the  resistivity  that  exists  across  the  width 
of  the  plasma  sheet  as  shown  in  Figure  II-2.  Note  that  the  positions  of 
minimum  resistivity  and  zero  magnetic  flux  density  along  the  y  axis  are 
not  identical.  Physically  we  can  see  the  mechanism  for  the  rippling  mode  by 
examining  the  three  regions  of  the  plasma  sheet  in  Figure  II- 2 .  Assume  first 
that  a  small  element  of  fluid  in  Region  l  is  given  an  upward  velocity  V  . 

This  carries  low  resistivity  fluid  into  a  high  resistivity  region.  Consequently 
a  first  order  current  density  will  flow  into  the  plane  of  the  drawing.  A 
motor  force,  J_j  x  B  ,  is  generated  that  opposes  the  upward  velocity,  V  . 

It  follows  that  this  region  is  stable.  Likewise  examination  of  region  3  shows 
that  it  is  stable.  Pushing  a  fluid  element  down  in  this  region  produces  a 
into  the  plane  of  the  drawing.  The  subsequent  x  B  is  upward.  In  Region  2, 
however,  an  element  of  fluid  moving  upward  carries  a  high  resistivity  into 
a  low  resistivity  region  creating  a  that  is  out  of  the  plane  of  the  drawing. 
The  motor  force  is  now  found  to  be  upward  along  the  y  axis.  This  region 
i.s  unstable. 
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Mathematically,  we  start  with  equation  2.  2)  and  expand  r|  and  J 
into  zero  and  first  order  components. 


r|  J  +  T)  J  =  V  x  B 
o—l  l—o  —  — 


2.  9) 


The  driving  force  for  the  rippling  instability  is 


F.  =  J  x  B  = - 

•—rippling  1 


Th  J  x  B 
1  — o  — 


V  B 


Since  the  above  equation  is  valid  for  a  small  region  ea  about  B  =  0  ,  the 
second  term  on  the  right  is  dropped.  We  now  have 


-rippling 


T)  J  xB 

I  — o  — 


2.  10) 


The  resistivity  r|  is  now  assumed  to  obey  the  continuity  equation 


below . 


5r]  +  V  ■  Vn  =  0 
M 


Expanding  and  retaining  only  the  first  order  terms  we  have 

tun,  +  v  V  =  o  2.ii) 

1  y  o 

Solving  for  and  rewritting  equation  2.  10)  we  get  the  following  result  . 
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But 

so  we  have 
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The  driving  power  for  the  mode  i6 


P-  V-F  .  =  -X-4  B'B 

-  -  rippling  M0'i,Tl0 
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This  power  must  be  of  the  same  order  as  the  power  given  in 
equation  2.  5)  in  order  to  generate  an  instability  at  all.  Therefore 


v2  n  '  2 

(Ca)  = 

U  't’Tl 
*o  o 


V2  (B')2  (ea)2 


2.  14) 


Equating  equations  2.  8)  and  2.  14)  we  get  an  expression  for  r  ,  which  is 
the  growth  rate  of  the  instability. 
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2.  15) 


The  derivations  of  F.K.R.  for  resistive  instabilities  employ  two 
characteristic  times  to  normalize  the  equation*.  They  are  : 


uo  a 

t  =  resistive  diffusion  time  «  - 

R  r\ 

o 

t  =  magnetohydrodynamic  transit  time  *  a 
H 


(S0‘ 


Equation  2.  15)  can  be  rewritten  in  terms  os  these  characteristic  times 


as  follows 


2  tr,  '\V5 

■ird 


u  (ka)2  (bVB)2  11/5 
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2  2  2 
47T  t,  a  u  C/B  J 


2.  II.) 


combining  terms 


whero  the  remainder  of  the  normalized  quantities  is  assumed  to  be  in  the 

order  of  unity.  This  result  is  in  agreement  with  the  detailed  calculations 

of  F,  K.  R.  For  convenience  we  finally  w rite  the  growth  rate,  normalized 

to  t0  ,  in  the  following  form. 

R 


where 
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2.  17) 


a 


The  "tearing"  mode  differs  from  the  rippling  mode  principally 
because  the  driving  force  for  tearing  is  not  concentrated  in  a  region  ea 
about  B  =  0,  but  depends  on  the  magnetic  field  structure  outside  the  region. 
The  mechanism  for  tearing  can  be  understood  from  an  examination  of 
Figure  II- 3  and  the  discussion  that  follows.  If  an  element  of  fluid  is  push(  1 
downwards,  the  magnetic  field  is  perturbed  as  shown  resulting  in  a  time  rate 
of  change  for  the  first  order  magnetic  field.  Taking  the  y  component  of  the 
following  equation, 


VxE 


SB 

=  '  TT 


we  have  -jk  E  y  =  j  B  y 

z  ■4-o  y  ■‘■o 

Thus  a  first  order  electric  field  is 

(i'E 


generated  that  has  the  following  form 
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The  interaction  of  this  electric  field  with  the  main  mangetic  flux 
density  produces  a  drift  velocity  downwards  in  the  direction  of  the  initial 
perturbation.  This  situation  is  unstable,  so  a  tearing  mode  develops. 


V  = 


E  x  B 


B 


uB 

- y 

jk  B  yo 
J  o  — 


The  electric  field  given  by  equation  2.  18)  is  complemented  by  a  first 
order  current  density,  J^.  Now  Ohm's  law  has  the  following  character. 

r  B 

2.19) 


n  J=E  ’•  V  x  B  =  -  — Z  +  V  x  B 
o—l  —  —  jk  _o  —  — 


For  the  instability  to  develop  the  first  term  on  the  right  must  dominate 
the  second  term. 
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it  B^ 

jk 


2.  20) 


The  value  Jj  can  be  related  to  the  magnetic  field  by 


V  x  B  =  u  J 

—  o—l 


and 

we  find  that  J 


V  •  B  =0 


B" 


and  if  B  is  approximated  bv  the  following 

juok  y 
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Then  we  find  from  substitution  in  equation  2.20)  that 
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_ 2 _ 1 _ 
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and  solving  for  ca  we  have 
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Equating  equations  2.8)  and  2.22)  an  expression  for  growth  rate 
of  the  tearing  mode  can  be  written. 


"1  ,4/5  ,2  ,2  1/5 
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Normalizing  via  the  characteristic  times  as  done  before,  the  above  equation 
becomes, 
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for  the  tearing  mode  where  p,  0.,  and  S  are  given  by  the  list  under  equation 
2.  17)  . 

The  second  part  of  this  chapter  will  be  concerned  with  the  development 
cf  theoretical  growth  rate  for  a  Rayleigh- Taylor  instability  along  a  plasma- 
vacuum  interface.  The  method  of  analysis  follows  that  c  ^senbluth  and 
Longmire  . 

We  shall  assume  an  infinite  plasma  carrying  current  with  a  magnetic 
field  supporting  the  plasma  against  gravity.  (Note  that  in  our  case  the 
acceleration  of  the  plasma  under  the  action  of  the  magnetic  pressure 
generates  an  equivalent  gravitational  force  .  )  It  will  be  shown  that  if  *he 
interface  is  perturbed  slightly,  an  instability  occurs.  Figure  II-4  shows  a 
sketch  of  the  situation. 

The  gravity  produces  a  drift  velocity  of  the  following  form 

V_  =  m  G  x  B 

q  b 


2.  24 


u 


so  that  the  ions  and  electrons  drift  in  opposite  direction?.  When  the  inter¬ 
lace  has  a  wavelike  perturbation  this  drift  results  in  charge  separation. 

This  separation  produces  an  electric  field  as  shown  in  the  figure.  This 
electric  field  oroducos  a  new  drift  velocity. 

V  -  i-Li  2.2 

£  B 

This  velocity,  V  ,  is  in  the  direction  that  tends  to  make  the  perturbation 
grow.  In  order  to  evaluate  the  growth  cate  we  assume  that  the  perutrbation 
is  of  the  form 


L  y  =  asm  kx 


2.  26) 


so  tiiat  the  change  in  charge  density, 6  ~,is 


6  r-  =  nq  6 Ay  =  nq  ^  fx 


t  ~  =  nq  a  k  cos  kx  fx 


2.27  ) 


,  then  we  have  for  fx  the  following 


m.  G 

vx  =  V  ft  =  ft 

F  qB 


2.  28) 


Here  the  electron  drift  nas  been  neglected  since  the  ion  mass  is  much 
greater  than  the  electron  mass. 

The  time  rate  of  change  for  the  surface  charge  density  -  (y)  is 

t  -  m  .G 

7 —  =  nqak  - - —  cos  kx 

ft  M  q  B 


6c  _  _i _ 

ft  B 


.  a  k  G 


cos  kx 
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where  =  nm. 

l  l 
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F^r  small  perturbations  and  a  very  thin  surface  charge  layer,  c  may  be 
expressed  as 

"  =  cr  cos  kx  2.  30) 

o 

The  electric  field  produces  by  this  charge  density  in  the  static  limit  can  be 
found  by  solving 

V  •  E  =  0 

v  x  E  =  -  —  =  0 
—  5,t 


The  electric  fields  in  the  plasma  are  thus 
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The  drift  velocities  associated  with  the  electric  fields  given  above  are 

E  "  lr 
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The  V  velocity  gives  rise  to  instability  and  at  the  boundary  y  =  0  causes 
the  perturbation  £y  to  grow  as  follows 
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Differ*  ntiating  equation  2 .  32)  again  we  get 
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so  that 

dfa 

2 

dt 

?.a  kG 
i 


B 


from  equations  2.  29)  and  2.  30) 
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From  equation  2.  33)  we  see  that  for  the  case  of  an  instability,  a  has 
the  following  form.  ^ 

a  =  exp  ( - -  kG  ^  t  2.34) 
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So  equation  2.  34)  becomes 

a  =  exp  v  kGt  2.35) 

where  G  does  not  necessarily  have  to  be  gravitational  acceleration  but 
can  be  a  result  of  magnetic  pressure  or  curvature  of  a  current  carrying 
plasma  . 

This  theory  neglects  any  damping  mechanisms.  One  form  of  dampi  lg 
can  be  a  result  of  viscous  Iriction  resulting  from  motion  of  the  plasma  under¬ 
going  a  Rayleigh- Taylor  instability.  The  damping  due  to  finite  viscosity  can 
be  calculated  simply  from  the  following  equation. 


{  •  (.-  •  V) 
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3<>) 


where  v  is  the  viscosity  . 
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The  plasma  configuration  is  the  same  as  shown  in  Figure  II- 4. 
We  assume  that  near  the  interface  the  velocity  has  the  forn 


-mt  Jky 


ko  . 


Substituting  in  equation  2.  36)  we  have 

mV  =  v  k2  V  2.  37) 

2 

Therefore  the  viscous  damping  constant  is  v  k  . 

Therefore  as  the  wavelength  of  the  instability  decreases,  damping 
due  to  finite  viscosity  increases  so  that  there  is  a  minimum  wavelength 
at  which  tile  instability  can  grow.  The  approximate  wavenumber,  k, 
correspond  ig  to  this  situation  can  be  determined  from  the  following 
equations, 


Solving  for  k,  we  have  : 


,2  „  J/2 

vk  =  (kG) 


,  _  1/3  -2/3 

k  =  G  v 


and  the  minimum  wavelength  is 


,  .  2tt  _  - 1/3  2/3 
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2.  38) 


2.  39) 


2.40) 


We  see  that  although  G  is  responsible  for  the  Rayleigh- Taylor  instability, 
its  influence  on  Amin  is  much  smaller  than  that  of  v,  and  therefore  >min 
is  essentially  a  measure  of  the  viscosity. 
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Development  of  Rayleigh- Taylor 
Instability 


Figure  Ii-4 
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CHAPTER  III 

ARRANGEMENT  OF  EXPERIMENTAL,  TOROID 

The  experiments  were  performed  on  an  electrodeless  toroidal  plasn  a 
oi  a  type  conceived  by  K.  T.  Eian^  for  the  purpose  of  high  power  switching. 

The  container  for  the  plasma  was  a  pyrex  toroid  having  a  rectangular  cross 
section.  During  the  course  of  the  experiment  both  narrow  and  wide  bottles 
were  used.  Shown  in  Figure  III  -  1  is  a  sketch  of  the  wide  toroidal  bottle. 

The  narrow  bottle  also  had  a  small  glass  pipe  inserted  radially  and  equi¬ 
distant  from  the  side  walls  of  the  toroid  so  that  magnetic  field  measurements 
in  the  core  of  the  discharge  region  could  be  obtained. 

The  toroidal  bottle  was  connected  to  a  high  vacuum  system  and  baked 

out  for  twenty  four  hours  at  400°  centigrade.  Vacuum  pressures  cf  about 

-8 

10  mm.  of  Hg  were  achieved  routinely  between  experimental  runs.  Hydrogen 
was  used  for  the  experiment  mainly  and  filling  pressures  varied  between  50 
microns  and  500  microns  of  Hg  depending  on  whether  Rayleigh- Taylor  or 
resistive  instabilities  were  being  studied. 

A  plasma  discharge  was  induced  in  the  hydrogen  by  means  of  an  iron 
cored  transformer.  A  hypersil  core  threaded  the  bore  of  the  toroidal  bottle 
and  a  primary  winding  was  placed  between  the  core  and  the  inner  surface,  of 
the  glass  envelope.  The  primary  winding  was  energised  by  a  high  voltage 
modulator.  When  the  hydrogen  in  the  bottle  becomes  ionized,  it  becomes 
a  gaseous  conductor  that  assumes  the  role  of  a  secondary  winding  and  load 
on  a  transformer. 

F’igui  e  III  -  Z  shows  the  various  primary  winding  configurations  used 
during  the  course  of  the  experiment.  The  function  of  these  different  winding 
arrangements  is  to  introduce  a  spatial  periodicity  in  the  electromagnetic 
field.  This  causes  an  initial  perturbation  in  the  plasma  that  determines  the 
wavelength  of  the  resistive  instability  that  will  develop.  Copper  strips  are 


23 


placed  over  the  windings  along  the  major  axis  so  that  perturbations  in  the 
applied  electric  field  are  minimized,  while  the  major  perturbation  is  due  to 
the  magnetic  field  developed  by  the  various  windings.  The  uniform  winding 
gives  rise  to  a  plasma  structure  with  two  filaments;  in  general,  if  N  is  the 
number  of  groups  into  which  the  primary  winding  is  divided,  the  discharge 
will  break  into  N  f  1  separate  filaments.  The  windings  are  fabricated  so 
that  the  ampere-turns  of  each  segment  are  cumulative  while  at  the  same 
time  generating  a  spatial  periodicity  that  causes  the  plasma  to  form  the 
appropriate  number  of  filaments. 

The  high  voltage  pulse  modulator  is  capable  of  delivering  up  to  20 
megawatts  power  to  the  toroidal  plasma  configuration.  Figure  III- 3  shows 
a  simplified  schematic  of  the  modulator.  Initially  the  R-L-C  combination 
is  impressed  across  the  D.C.  voltage  power  supply.  The  L-C  elements 
form  a  tank  circuit  and  therefore  the  voltage  and  current  are  of  an  oscillatory 
nature.  Because  of  the  diode  shown  in  the  circuit,  only  the  first  h:<lf  cycle 
of  oscillation  is  allowed  and  the  voltage  drop  across  the  capacitor  at  the  end 
of  charging  is  double  that  of  the  D.C.  voltage  supply.  The  resistor,  R^, 
serves  to  limit  the  charging  current  while  R^  is  adjusted  so  that  the  pulse 
current  to  the  toroid  is  essentially  critically  damped. 

The  current  pulse  is  applied  to  the  toroid  by  triggering  the  hydrogen 
thyratron.  A  low  energy  pulse  with  a  60  Hertz  repetition  rate  is  applied  to 
a  driver  unit  that  amplifies  the  trigger  and  then  drives  the  grid  of  the  hydrogen 
thyratron  positive,  thus  firing  it  into  a  conducting  state.  The  current  pulse 
has  a  duration  of  approximately  six  microseconds  and  its  peak  amplitude 
was  adjusted  between  4  kiloamperes  and  7  kiloamperes  during  the  course  of 
the  experiment. 

Preionization  of  the  toroidal  plasma  was  employed  to  insure  that  the 
plasma  was  initially  uniform  in  the  axial  direction  before  the  large  current 
pulse  from  the  modulator  was  applied.  Energy  for  preionization  is  coupled 
into  the  plasma  capacitvely  via  the  two  copper  sheets  encircling  the  outer 
circumK  rence  of  the  toroid  as  shown  in  Figure  III- 1 .  An  overall  schematic 
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of  the  experiment  is  given  in  Figure  III-4.  In  this  figure  is  shown  the  R.  F. 
transmitter,  which  is  tuned  to  7  MHz,  that  supplies  power  to  a  tank  circuit. 
The  inductance  of  the  tank  circuit  is  about  1  20  u  h  wound  on  a  six  inch  coil 
form  and  this  is  connected  in  parallel  to  the  copper  sheets  that  form  the 
capacitor  mentioned  above. 

Also  shown  in  the  figure  is  the  modulator  that  suDplies  the  current 
pulse  to  the  toroid,  timing  circuitry,  and  various  instrumentation.  Included 
in  instrumentation  are  an  S.  T.  L.  image  convertor  camera  for  taking  ultra 
high  speed  photographs,  monochromater  and  spectrograph  for  making  optical 
measurements  of  plasma  electron  temperature  and  density,  and  probes  for 
measuring  magnetic  field,  plasma  current,  and  modulating  voltage.  Elec¬ 
trical  outputs  from  the  instrumentation  are  monitored  on  a  Tektronix  555 
dual-beam  oscilloscope.  Electrical  timing  of  instrumentation,  oscilloscope 
and  high  voltage  pulse  modulator  are  accomplished  by  a  60  Hz  clock  dri\  ng 
an  S.  T.  L.  pulse  generator.  This  generator  develops  fixed  and  time  delayed 
pulses  that  have  rise  times  of  2  nanoseconds.  The  time  delayed  pulses  can 
be  varied  between  0  -  100  microsecond  delay  time  in  10  nanosecond  steps 
with  an  accuracy  of  one  nanosecond.  This  results  in  an  ability  to  record 
data  with  great  precision  and  makes  correlation  of  various  data  from  the 
instrumentation  easier. 
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Primary  Winding  Configurations 


Figure  III- 2 
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CHAPTER  IV 


INSTRUMENTATION 


Several  methods  are  employed  to  determine  the  properties  that 
characterize  the  toroidal  plasma.  Voltage  and  current  measurements, 
magnetic  probes,  fast  photography  with  an  image  converter  camera,  and 
spectroscopic  determination  of  density  and  temperature  are  included  in 
the  instrumentation  methods. 

A.  Voltage  Measurements 

Voltage  measurements  are  obtained  with  a  Tektronix  P6015  high  voltage 
probe.  The  prcbe  is  attached  to  a  single  turn  of  wire  encircling  the  magnetic 
core  of  the  toroidal  plasma  bottle.  The  voltage  measured,  then,  is  equivalent 
to  the  voltage  drop  for  a  complete  circular  path  ii.  the  plasma  ring. 

The  voltage  signal  is  helpful  for  determining  rapid  inductive  changes 
in  the  system.  For  example,  at  breakdown  of  the  gas,  the  initially  high 
voltage  drops  to  a  low  value.  Also  changes  'n  radial  plasma  motion  can  be 
detected  by  the  behaviour  of  the  voltage  signal. 

B.  Current  Measurements 

The  current  fed  to  the  primary  winding  from  the  modulator  is  monitored 
by  a  current  transformer.  The  current  transformer  consists  of  a  permalloy 
core  about  1  1/2"  in  diameter  and  the  primary  of  this  transformer  is  a  single 
wire  passing  through  the  center  of  the  core.  This  wire  carries  the  modulator 
current.  The  secondary  winding  consists  of  many  turns  coupled  to  a  resistive 
load  of  approximately  a  few'  ohms.  The  voltage  developed  across  this  resistor 
is  proportional  to  the  modulator  current  if  the  leakage  flux  of  the  transformer 
is  neglected. 

Actually  it  is  the  measurement  of  plasma  current  that  is  of  prime 
interest.  Although  this  current  can't  be  measured  directly,  the  difference 
between  primary  current  and  plasma  current  can  be.  A  core  for  a  difference 
current  transformer  is  made  of  permalloy  tape.  This  tape  is  threaded  so 
that  it  links  both  the  primary  winding  and  the  plasma  in  the  toroidal  bottle  . 
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The  ampere-turns  associated  .with  the  primary  winding  and  the  plasma  make 

up  the  primary  flux  for  the  difierence  current  transformer.  Again  many 

turns  encircling  the  tape  form  a  secondary  winding  and  this  is  terminated 

in  a  resistor  of  only  a  few  ohms.  Ideally  the  voltage  developed  across  this 

resistor  should  be  zero;  for  this  indicates  perfect  coupling  between  the 

plasma  and  the  primary  winding  of  the  toroidal  bottle.  However,  if  the 

voltage  is  not  zero  the  plasma  current  can  be  computed  since  the  difference 

current  transformer  can  be  calibrated  against  the  modulator  current  in  the 

absence  of  a  plasma,  and  then  the  plasma  current  is  computed  as  the 

following:  I  =N.  I  — I  , 

plasma  primary  primary  dm.  cur.  xfmr. 

C.  Magnetic  Probes 

Both  stationary  and  movable  probes  were  constructed.  The  stationary 
probes  were  mounted  on  two  lucite  disks.  Each  disk  was  positioned  on  a  side 
wall  of  the  toroidal  bottle.  Imbedded  on  each  disk  are  ten  concentric  loops 
as  shown  in  Figure  IV- 1.  Their  main  purpose  is  to  sample  the  axial  com¬ 
ponent  of  the  magnetic  flux  density,  Bz,  and  thus  determine  its  time  re¬ 
solved  variation  of  the  radius.  From  this  variation  one  can  ascertain  the 
radial  motion  of  the  plasma  current  sheet.  The  symmetry  and  stability  of 
this  motion  is  checked  by  four  small  probes  mounted  on  each  lucite  disk 
90°  apart.  These  probes  sample  Bz  also.  Further  examination  of 
Figure  IV- 1  reveals  another  set  of  small  fixed  probes  that  are  orthogonal 
to  the  set  of  probes  mentioned  above.  This  last  group  of  four  probes  are 
also  90°  apart  and  measure  the  symmetry  of  the  radial  component  of  magnetic 
flux  density,  Br. 

The  small  movable  probe  is  shown  in  Figure  IV-2.  Two  such  probes 
were  made,  one  to  measure  Br,  and  the  second  to  measure  Bz  (B«-can  be 
measured  by  merely  rotating  this  last  probe  90°).  The  probes  consist  of 
a  coil  of  50  turns  of  #50  enameled  copper  wire  wound  on  a  ~  inch  bakelite 
rod.  These  probes  can  be  moved  around  and  inside  the  glass  envelope  of 
the  plasma  toroidal  bottle  allowing  the  magnetic  field  to  be  mapped.  Appli¬ 
cation  of  Maxwell's  equations  yields  the  current  distribution. 
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The  leads  from  all  probes  were  twisted  to  eliminate  stray  pickup 
In  addition  the  small  movable  probes  incorporated  both  magnetic  a"d 
electric  shielding  on  the  leads.  Signals  from  the  probes  were  electronically 
integrated  v,  ith  a  RC  circuit  to  obtain  Br,  Bz,  and  B-e-.  The  time  constant  of 
tlie  RC  circuit  was  80  usee.  ,  which  is  approximately  ten  times  the  modulator 
pulse  length. 

D.  Fast  Photography  with  a  S.  T.  L.  Image  Converter  Camera 

An  independent  sampling  of  the  dynamics  of  the  plasma  discharge  is 
obtained  with  the  help  of  a  S.  T.  L.  image  converter  camera.  This  camera 
is  fitted  with  two  tubes,  one  with  sensitivity  peaking  in  the  red  and  the  other 
peaking  in  the  blue  portion  of  the  optical  spectrum.  Both  tubes  have  a  light 
flux  gain  greater  than  50  X  each.  It  is  thus  possible  to  obtain  a  photograph  of 
the  overall  light  emission  of  the  plasma,  and  by  using  appropriate  inter¬ 
ference  filters,  to  follow'  the  space  and  temporal  development  of  specific- 
spectral  lines. 

The  camera  can  be  operated  in  two  modes:  framing  and  streaking  . 

In  the  framing  mode  a  sequence  of  nine  pictures  of  the  ^asma  discharge 
can  be  taken  with  an  exposure  time  of  50  nanoseconds  typically  and  with  a 
time  interval  between  pictures  ol  500  nanoseconds.  When  operating  in  the 
streaking  mode,  the  camera  is  u^e.  '  ..  a  lusec/  10mm.  writing  speed. 

However  the  camera  can  scan  as  slow  as  1  usec/5mm.  or  as  fast  as  1  nanosec- 
/lmrn.  The  camera  is  triggered  by  a  variable  time  delay  pulse  from  a  S.  T.  L. 
trigger  delay  generator.  The  same  generator  triggers  the  thyratron  in  the 
modulator  with  a  fixed  pulse.  Satisfactory  time  resolution  in  the  photograph 
sequence  is  assured,  cince  the  jitter  in  the  time  delay  doesn't  exceed  Z  nano¬ 
seconds  . 

Figures  IV- 3  and  IV-4  show  the  camera  set  up  to  take  a  sequence  of 
framing  pictures  along  the  z-axis  and  perpendicular  to  the  z-axis  respectively. 
When  taking  streak  photographs  it  is  necessary  to  select  light  emission 
whose  rays  are  almost  parallel,  in  other  words,  a  small  solid  angle. 
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Streak  photographs  taken  of  light  along  the  radial  axis  require  no  auxiliary 
equipment  since  the  camera  must  be  far  from  the  plasma  bottle  in  order  to 
accept  all  radial  light  rays  emitted  along  the  z  axis  of  the  plasma  discharge. 
See  Figure  IV--5  .  However  streak  photographs  of  2  axis  light  emission  with 
a  radial  distribution  have  a  smaller  field  of  view  than  in  the  previous  case. 
This  requires  the  camera  to  be  closer  to  the  event  which  causes  the  solid 
angle  to  increase.  In  this  case  a  mechanical  collimator  with  a  1/2°  accept¬ 
ance  angle  Is  mounted  directly  in  front  of  the  camera.  This  set-up  is  illus¬ 
trated  in  Figure  IV-6  . 

E.  Quantitative  Streak  Photography 

Streak  photographs  are  particulary  suited  for  quantitative  measurements 
of  variation  in  light  intensity  emitted  by  the  plasma  because  this  photograph 
is  a  continuous  temporial  record  of  either  radial  or  axial  distribution  of 
essentially  parallel  light  rays.  A  major  difficulty  is  the  fact  that  the  infor¬ 
mation  is  stored  logarithmically  and  a  suitable  calibration  of  the  film  must 
be  ni»Je. 

The  method  employed  here  is  to  reserve  a  small  portion  of  the  film 
for  making  a  calibration  record.  The  plasma  bottle  is  completely  masked 
except  for  a  small  pinhole.  A  ccreak  photograph  of  light  emission  from 
pinhole  is  recorded  in  the  reserved  spot.  Next  a  photodiode,  which  has  a 
linear  response  to  light  intensity,  is  exposed  to  the  pinhole  .  The  photo¬ 
diode  signal  is  fed  to  an  oscilloscope  and  photographed. 

A  Wihtol  microdensitometer  and  Brown  Electronik  recorder  shown 
in  Figure  IV-7  are  used  to  analyse  the  film.  First,  the  calibration  record 
on  the  film  is  scanned  and  plotted  against  the  photodiode  signal  at  identical 
instants  of  time.  This  results  in  a  calibration  curve  for  the  film.  Then  the 
streak  pictures  of  the  plasma  are  scanned,  each  data  point  being  unfolded 
by  the  calibration  curve  to  its  proper  value  of  light  intensity. 

This  technique  was  used  extensively  for  the  determination  of  growth 
rates  for  resistive  instabilities  . 
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F .  Spectrograph 

A  set-up  employing  a  Baush  and  Lomb  1.  5  meter  spectrograph  to 
survey  the  spectrum  emitted  by  the  plasma  discharge  is  shown  in  Figure 
IV-8.  Royal  X  pan  film  was  used.  An  exposure  time  of  approximately 
15  minutes  and  development  in  Acufine  was  required  to  produce  acceptable 
results.  The  film  was  then  microdensitometered  in  order  to  facilitate 
interpretation  of  the  spectrum  for  impurities.  In  addition  the  Balmer 
series  spectral  lines  or  Hydrogen  were  examined  and  the  number  density 
calculated  using  the  method  of  Inglis  and  Teller.  In  this  technique  the  micro 
fields  of  both  ions  and  electrons  cause  broadening  of  the  higher  energy  levels 
in  the  series.  When  the  last  distinguishable  line  has  a  width  half  the  "distance" 
to  the  preceding  line  the  relation  between  electron  density  and  the  last  line 
is  the  following 
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where 


a  -  Bohr  radius 
o 

Z  -  atomic  number 

n  -  total  quantum  number  of  last  distinguishable  line 
m 

e  -  electron  charge 
N  -  electron  density 


The  electron  density  can  thus  be  ascertained  by  identifying  the  last 
distinguishable  line. 
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G.  Monochromate  r 

A  500mm.  Baush  and  Lomb  grating  monochromate r  equipped  with 

a  photomultiplier  tube  was  used  to  examine  the  iight  emission  in  detail. 

o 

The  monochromater  was  attached  to  a  motor  drive  so  that  scans  of  2A  / 
o 

second  or  20  A /  second  could  be  made.  In  addition, Hilger- Watts  slits 

(repeatability  of  less  than  2  microns)  were  used  to  increase  the  systems 

efficiency.  Following  Griem, ^  temperature  measurements  were  made 

by  comparing  the  ratio  of  the  H.  line  in  the  Balmer  Series  of  Hydrogen  to 

o  0 

the  neighboring  100  A  continuum.  This  measurement  was  performed  at  a 
o 

2A /second  grating  speed,  and  since  the  plasma  discharge  is  pulsed,  a 
Tektronix  151  plug-in  was  used  to  sample  the  photomultiplier  tube  output 
at  a  particular  instant  of  time.  The  output  of  the  151  plug-in  unit  drove  a 
strip  chart  recorder  (see  Figure  IV-9)  .  Thus,  the  plot  on  the  chart  was  a 
time  resolved  spectrogram.  This  set-up  was  prepared  by  R.  Chimenti. 

o 

A  planimeter  was  used  tc  measure  the  area  under  the  H.  line  and  100A 

6  o 

continuum  so  that  the  ratio  could  be  calculated.  The  curve  of  H./100A 

o 

continuum  versus  temperature,  taken  from  Griem,  is  shown  in  Figure  IV- 10. 

A  technique  for  measuring  electron  number  density  via  examination 

12 

of  Stark  profiles  described  in  a  paper  by  Griem  et.  al.  appeared  attractive. 
In  this  method,  the  Stark  profile  is  not  measured  per  se,  instead  the  mono¬ 
chromater  is  set  at  the  line  core,  and  photomultiplier  tube  outputs  are 
recorded  for  narrow  and  wide  slits.  The  ratio  of  narrow  to  wide  exit  slit 
signals  uniquely  determines  the  number  density  and  this  measurement  is 
practically  independent  of  temperature.  Griem  worked  with  H.  ;  however 
in  my  experiment  H  was  selected  as  the  optimum  choice.  The  curve  of 

U 

density  versus  the  ratio  of  signals  with  ,  1mm  exit  slit  and  1mm  exit  slit 
is  shown  in  Figure  IV- 11  . 
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Tht  agreement  between  the  experimental  line  profiles  and  those  due 
to  the  general  theory  of  line  broadening  set  forth  by  Griem,  Kolb  and  Shen 
is  about  5%  for  II  .  Unfortunately  the  solutions  for  the  theoretical  line 
profiles  can  not  be  expressed  in  closed  form,  so  numerical  tables  have 
been  prepared  by  the  authors  which  enable  one  to  plot  a  particular  line 
for  various  electron  densities  and  temperatures. 

14 

The  Holtsmark  theory,  presented  in  1919,  yields  a  profile  that 
gives  good  agreement  in  the  wings,  and  moreover  predicts  the  correct 
half  width  of  the  broadened  line. 


The  sketch  shown  above  illustrates  the  comparison  between 

1  3 

Holtsmark's  theory  and  that  due  to  Griem  et.  al.  Disagreement  between 
the  two  curves  is  pronounced  in  the  neighborhood  of  the  line  core.  He. 
doppler  broadening  plays  a  greater  role.  This  phenomena  is  nof  taken  into 
account  by  the  Holtsmark  theory. 

o 

The  following  is  an  explanation  of  why  the  ratio  of  light  flux  for  1.  bA 
o 

and  16A  wavelength  bands  centered  at  the  line  core  of  the  Hq  line  results 
in  a  unique  determination  of  electron  density.  First  it  is  asserted  that  the 
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o 

effect  of  dooplcr  broadening  is  confined  mostly  in  the  1. 6A  wavelength 

r> 

band,  and  for  larger  wavelength  ba'.ds,  eg.  16  A  the  Holtsrr  irk  theory  is 
a  fair  approximation  to  the  actual  line  profile,  S(AA)  .  We  next  assume  that 
,he  total  line  prof’lc  is  prescribed  by  the  Holtsmark  function  given  below' 
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A^,  shift  in  wavelength  of  k —  line  due  to  Stark  effect  . 
k 

=  constant  dependinb  on  initial  and  final  quantum  numbers 
involved  in  the  transition  for  excitation  of  H  . 

D 

E  =  "average"  electric  field  associated  with  the  Stark 
broadeninp  effect  . 


The  light  intensity  distribution  for  (he  Stark  broadened  line  is  related  to  the 


Holtsmark  function  by 

£ 

Tk  =  c  H<B> 


4.  4) 


where 


1^  =  intens  ty  of  unbroaderad  k—  line 
1^  =  intensity  of  Stark  broadened  k~  line 


Now  equations  4.  3)  and  4.  4)  can  be  combined  to  yield 
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The  average  electric  f'eld,  E  ,  is  r  ’ally  a  function  oi  tha  electron 

o 

density.  This  can  be  demonstrated  by  first  calculating  the  verage" 
separation  between  ions  for  a  density  of  ions  equal  to  N.  .  7  his  is  done  by 

solving  the  equation  be  ow. 
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where  R  =  average  reparation  ,  N.  =  ion  density  . 
electric  field  is 


Now  the  average 


E  = 
o 


4tt  e  R2 
o  o 


4.  7) 


Substituting  equation  4.  (>)  into  4.  7)  we  have 


E  =  AN  2A 
o  1 

where  A  =  constant 

For  hydrogen  we  have  electron  density  equal  to  icn  density,  therefore 

E  =  AN  % 
o  e 

where  N  =  electron  density, 
e 

finally  the  substitution  of  equation  4.  8)  into  4.  5)  results  in 
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where  C.  =  AC, 
k  k 


The  total  light  flux  of  the  Stark  broadened  spectral  line,  over  a 

wavelength  band  AX  is  the  following: 
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If  we  now  let  A X  ^  =  l.bA  wavelength  band  and  -  1  6A  wavelength 
band  and  form  the  ratio  of  light  flux  for  the  two  bands 


have  : 
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Simplifying. equation  4.  1 1 )  becomes 
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Nov;  the  variation  of  the  light  flux  ratio,  R,  with  electron  density  is 
given  by 


dR  _  —  i  »•  |  . . 

dN  '  "  °k  VAX1 
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AX2  -  J  R 


4.  13) 


We  know  that  C.  >0 
k 

0<  R <  1 


and 


AX,  <AX, 


\<  ■%  -  A X  ~ %  >0 


The  above  implies  that  <0  ,  or  that  R  monotonically  decreases  as 

N  increases.  This,  in  fact,  is  the  behavior  of  Figure  IV-11  which  shows 
e 

N  versus  A,  /A,  . 
e  12 

Electron  density  was  also  determined  by  inserting  the  observed  half 
width  of  a  line  into  an  equation  of  the  form 


N  =  C  AX  % 
e 

where  C  is  a  quantum  mechanical  constant  that  is  a  weak  function  of 

15 

electron  density.  Such  an  expression  has  its  roots  from  the  equation 
given  below,  which  relates  the  wavelength  shift  of  a  spectra)  line  to  ar. 
impressed  e^ctric  field.  This  is  Ihe  first  order  Stark  effect. 


4.  14) 
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L  X  ~  X^  (nn,_  -  n'n  )  E  4.15) 

F  F 

where 
and 


n,  n7  =  p  ncipal  quantum  numbers 

n  ,  n^,  =  electric  field  quantum  numbers  related  to  electric 

dipole  moment  of  atom  and  oerturbing  electric  field,  E. 


This  last  equation  has  the  same  form  as  equation  4.  14)  .  Values  of  C  are 
given  in  the  book  "Plasma  Spectroscopy"  by  Griem. 
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OF  THE  SO  TURN 
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Stationary  Magnetic  Probes 
Figure  IV-  1 
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Spectrograph  Setup 
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CHAPTER  V 


EXPERIMENTAL  FINDINGS 


This  chapter  presents  the  results  that  were  obtained  by  the 
instrumentation  described  in  Chapter  IV.  First,  the  measurement  of 
magnetic  fields  is  given.  Next  the  experimental  results  are  presented  for 
measurement  of  electron  density  and  temperature.  The  following  section 
gives  the  values  of  resistivity  for  various  modulator  voltages  and  plasma 
currents.  The  last  two  sections  are  concerned  with  the  observation  and 
measurement  of  growth  rate  for  the  Rayleigh- Taylor  instability  and  the 
resistive  instability. 

A.  Measurement  of  Magnetic  Fields 

Detailed  measurements  of  the  magnetic  flux  density  were  carried 

out  on  the  narrow  toroidal  plasma  bottle.  The  first  measurements  made 

on  this  type  of  plasma  structure  was  done  by  H.  Friedman,  ^  who  compared 

his  data  to  a  theoretical  plasma  model  that  consisted  of  a  current  sheet. 

Lat  r,  a  more  detailed  study  of  the  magnetic  flux  dt  nsity  was  undertaken 
17 

by  J.  Milleta  .  He  also  compared  the  experimentally  measured  fields  <vith 
a  model  that  represented  the  current  carrying  plasma  as  being  comprised 
of  two  rings.  This  model  was  a  direct  result  of  the  examination  of  framing 
photographs  such  those  given  in  Figure  V-29  which  show  the  plasma 
structure  in  the  narrow  toroidal  bottle  being  comprised  of  two  rings  or 
filaments. 

The  magnetic  probes  used  were  Df  two  types,  stationary  and  me '/able. 
The  stationary  probes  were  mouw.ed  on  the  side  walls  of  the  plasma  bottle 
and  the  movable  probes  had  the  facility  of  measuring  magnetic  flux  density 
both  inside  and  outside  the  toroidal  bottle  .  The  probes  wer  arranged  so 
that  all  three  components  of  magnetic  field,  Br,  B-n-,  and  B  c  ould  be 
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measured.  The  complete  description  of  the  magnetic  probes  and  their 
position  has  been  given  in  Chapter  IV-C. 

Shown  in  Figures  V-l,  2,  3  are  photographs  of  temporal  variations  of 
Br,  Boi  Bz  respectively  for  different  radial  positions  equidistant  from  the 
side  walls  of  the  toroidal  bottle.  Also  shown  for  reference  purposes  is  the 
primary  current  as  a  function  of  time.  Detailed  magnetic  prebe  measure¬ 
ments  were  also  made  along  the  axial  length  of  the  outer  circumference  of 
the  toroidal  bottle. 

These  data  for  the  three  components  of  magnetic  flux  density  are 
functions  of  time  and  space.  To  aid  in  the  interpretation  of  this  large 
amount  of  data,  plots  cf  flux  lines  were  made  at  instants  of  time  corres¬ 
ponding  to  one,  two,  three,  and  four  microseconds  after  the  beginning  of  the 
current  pulse,  The  plots  were  made  for  a  cross-section  of  *-He  toroidal  bottle, 
the  r ■  z  pLrne.  Also  shown  in  the  plots,  in  out  line  form,  is  the  perimeter  of 
the  bottle,  and  also,  the  magnetic  flux  density  vectors  that  were  calculated 
from  the  experimental  data.  The  plots  depicting  lines  of  flux  for  various 
instants  of  time  are  shown  in  Figures  V-4  through  V-7  inclusive. 

Examination  of  these  figures  reveals  two  separate  flux  linkages  , 
which  implies  that  the  plasma  structure  consists  of  two  current  carrying 
rings.  Moreover,  the  two  flux,  patterns  shown  for  each  instant  of  time  are 
not  identical.  This  suggests  that  one  plasma  ring  conducts  more  current  than 
the  other.  This  observation  is  also  consistant  with  the  fact  that  one  current 
ring,  the  one  with  more  flux  linkages,  has  a  greater  outward  radial  motion. 
This  motion  can  be  related  to  the  well  known  "magnetic  piston"  argument 
which  stems  from  the  J[  x  B  motor  force.  If  the  current  in  the  ring  is 
larger,  then  the  Jx  B  force  increases, and  the  subsequent  outward  radial 
motion  is  faster. 
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Photographs  also  support  the  argument  of  separate  filaments.  In 
addition  f  Figure  V-29  mentioned  previously,  another  framing  photo¬ 
graph  that  shows  radial  motion  is  given  in  Figure  V-8  .  In  this  figure,  the 
camera  is  looking  along  the  z  axis  through  an  observation  window  in  a  side 
wall  of  the  toroidal  bottle,  The  time  sequence  is  given  in  the  figure.  Of 
particular  interest  are  the  photographs  for  instants  of  time  at  1.0,  1.5,  and 
2.0  microseconds.  Here  two  rings  can  be  seen,  one  of  them  being  at  a 
smaller  radial  distance  from  the  center  of  the  toroidal  bottle. 

The  current  sheet  model  and  the  two  ring  model  for  the  magnetic  fields 
in  the  toroidal  plasma  discharge  both  neglect  the  curvature  of  the  torodial 
structure  and  also  assume  that  the  radial  dimension  of  the  plasma  is  infini¬ 
tesimally  thin.  A  schematic  of  the  theoretical  models  is  given  in  Figure  V-9. 
Solutions  for  the  Br  and  Bz  components  of  magnetic  flux  density  using  a 
current  sheet  model  are  given  as: 

Br  =  [uoI{t)/2nd]  +  (z +d/2)2]/ fr^  +  (z-d/2)2]  | 

Bz  '  r-uoKt)/4TTd]  -jtan  *  [  (z  +  d/2)/r]  -  tan  *  ["  (z-d/2)/r  ]  j 

The  total  magnetic  field  is  the  superposition  of  the  fields  due  to  the  plasma 
current  sheet  and  the  primary  current  flus  their  image  fields,  where  it  has 
been  assumed  that  the  magnetic  transformer  core  is  a  semi-infinite  magnetic 
wall.  A  sketch  of  this  appears  in  the  bottom  portion  of  Figure  V-9. 

The  same  prescription  is  followed  for  the  calculation  of  total  magnetic 
flux  density  in  the  case  of  the  two  ring  model.  The  equations  for  Br  and  Bz 
components  are  given  below.  The  distance  c?/  represents  the  axial  separation 
of  the  two  rings.  It  is  assumed  that  both  rings  conduct  equal  current. 

Br  =  Tuo  I(t)/4rr]  |(z-dV2)/  [r^+  (z-dV2)l  +  (z  +  a  11)1  Tr^  +  (z  +  dV2)]  j  5.  3) 

Bz  =  f -Uq  I( t ) / 4tt 3  |  l/[r^  +  (z-d'/2)^]  +  1/  [r^  +  (z  +  d'/2)^]  j 


5.1) 

5.2) 


5.4) 


» 
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The  magnetic  fields  fur  both  models  are  plotted  for  the  following 
conditions  : 

modulator  voltage -  12kv 

neutral  gas  pressure -  500  microns  Hq 

plasma  current  -  3400  A;  magnitude  at  2  usee. 

after  beginning  of  current  pulse. 

Fir  ires  V-10  and  11  show  the  Br  and  Bz  components  due  to  the  models. 
Comparison  with  the  flux  plot  for  the  time  of  two  microseconds  indicate 
that  the  actual  field  i.*;  not  given  in  detail  by  either  of  the  models.  However, 
if  the  two  ring  model  consisted  of  unequal  currents  in  the  rings,  the  differ¬ 
ence  between  model  and  actual  field  distribution  would  be  less. 

Of  more  practical  interest  for  the  purposes  of  instability  calculations 
is  the  value  of  the  typical  magnetic  field  ,  <B>,  within  the  plasma.  It  was 
found  that  this  magnetic  flux  density,  <B>,  expressed  by 

4 

<B>  *  10  uo<I>/2L  (gauss)  where  5.5) 

L  is  the  axial  length  of  the  plasma,  matches  the  value  of  both  model  fields 
and  the  actual  field  associated  with  the  portion  of  plasma  in  the  center  of 
the  toroidal  bottle. 

The  experiments  on  instabilities  of  the  resistive  type  indicated  their 

presence  at  times  shortly  before  the  peak  of  the  current  pulse,  approximately 

two  microseconds.  At  this  time,  the  variation  of  magnetic  flux  density 

from  outer  plast.  a- vacuum  to  inner  plasma- vacuum  interface  was  110 

gauss  to  360  gauss.  This  represents  an  average  field  of  235  gauss.  For  an 

2 

L~  6  cm.  and  an<I>  =  —  I  peak,  where  I  peak  =  3400A  at  two  microseconds  , 
we  find  that  <B>  is  as  follows  : 


<B> 


<I> 

2L 


<B>  =  226  gauss 

Calculation  of  the  average  model  field  at  two  microseconds  yielded  a 
value  of  22G  gauss. 
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Therefore,  when  magnetic  flux  density  is  required  for  instability 
calculations,  the  approximation  expressed  by  equation  5.5)  v; ill  be  used. 


B.  Measurement  of  Electron  Densities 

Electron  densities  were  measured  in  the  tc-roidal  plasma  by  obser¬ 
ving  the  effect  of  Stark- broadening  on  the  Hydrogen  spectral  lines  in  the 

1 8 

Balmer  Series.  The  theory  of  Inglis  and  Teller  was  used  to  calculate 
den  -  from  the  experimental  observation  of  line  merging.  Several  line 
shapes  were  examined  in  detail,  and  comparison  of  these  shapes  with  the 
theoretical  ones  presented  by  Griem^  yield  the  electron  density. 

The  first  measurement  of  electron  density  employed  the  spectrograph 
set-up  that  is  discussed  in  Chapter  IV-F.  In  this  measurement  the  density 
was  not  time  resolved,  since  the  shutter  of  the  spectrograph  remained  open 
for  approximately  fifteen  minutes. 


A  portion  of  the  microdensitometered  film  appears  in  Figure  V-l.t. 
Line  merging  was  considered  to  be  significant  after  the  line  labeled  n=  10 
because  the  wings  of  all  succeeding  lines  overlap  to  a  large  extent.  Using 
the  formula  given  in  equation4.  2)  of  Chapter  IV  for  density  we  have 


N 


% 


1 


i,  .  5  2 

1 1 .  1  n  a 
m  o 


where  n  =  10,  and  a  =  .  53x  10  cm.  Solving,  we  have  N  =5.7x10 
m  o  B  e 

part/cm^  ,  This  data  was  obtained  with  the  modulator  set  at  10  kv  and 
with  a  neutral  gas  pressure  of  500  microns. 

Later,  a  time  resolved  scan  of  the  spectrum  was  achieved  using  the 
monochro: rater  set-up  as  discussed  in  Chapter  IV-G.  The  photomultiplier 
tube  output  of  the  monochromater  was  strobed  or  sampleu  at  4  microseconds 
from  the  beginning  of  the  current  pulse,  while  simultaneously,  the  mono¬ 
chromater  grating  was  being  rotated  by  the  motor  drive.  The  resultant 
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t  utput  is  shown  in  Figure  V-13  ,  this  figure  of  course,  is  only  a  portion  of 

o 

tne  total  spectrum.  The  scale  is  15.  6  A  /inc.i  -vid  only  the  lines  with 

principal  quantum  number  n=£  . Trough  n=10  are  shown.  Of  interest  is  that 

the  line  associated  with  n=10  is  absolutely  the  last  observable  line;  all 

jucceedi'ig  lines  are  buried  in  the  continuum.  For  n=10,  we  have  again 

15  3 

from  the  Inglis  and  Telle  '  technique  a  deraity  of  5.  7x10  part/cm  .  This 
value,  however,  is  time  resolved  at  4  microseconds. 

A  more  accurate  determination  of  number  density  is  to  measure  the 
half  width  of  a  Hydrogenic  spectral  line  or  to  plot  the  profile  itself,  these 
measurements  are  then  compared  to  the  theory  of  Griem.  The  mon*chro- 
mater  is  adjusted  so  that  it  is  set  at  the  line  core  of  H_  .  Then  the  photo- 
multiplier  tube  output  is  recorded  with  the  exit  slit  set  first  at  0  1mm  and 
the"  xt  1.  0mm.  Tl.  •  ratio  of  outputs  at  those  two  slit  openings  represents 
a  unique  electron  density.  Such  a  cu.ve  is  shown  in  Figure  IV-11. 

Data  was  taken  for  a  range  of  neutral  gas  pressures  ar.d  modulator 
voltages.  The  oscilloscope  photographs  or  photomulitplicr  tube  output  show 
signal  variation  with  time. 

It  was  noticed  that,  continuum  radiation  was  very  large  for  the  first 
microsecond,  with  its  peak  at  about  1/2  microsecond.  After  one  micro¬ 
second,  continuum  ’•adiation  was  negligible  when  compared  with  the  Ilg 
light  output.  For  t.mes  after  or.e  microsecond,  therefore,  continuum  is 
neglected  when  computing  density  from  the  profile.  For  shorter  times  , 

P 

the  continuum  radiation  must  be  subtracted  from  the  H0  data. 

s 

Figure  V-14  shows  the  following  data:  in  photograph  a  are  the 
primary  current  and  core  voltage  a  a  function  of  time;  photograph  b  shows 
photomultiplier  tube  output  with  narrow  slit,  and  photograph  c  is  taken  with 
a  wide- slit.  This  figure  represents  data  from  a  torodial  plasma  driven  at 
a  modulator  voltage  of  12kv  and  neutri'l  gas  pr«  ssure  of  500  nr  .tops. 
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Graphs  depicting  the  electron  density  as  a  function  of  modulator  voltage 
and  neutral  gas  pressure  for  instants  of  time  corresponding  to  two  '  <id  four 
microseconds  after  the  beginning  of  the  current  pulse  are  shown  in  Figure 
V-15  to  18  inclusive. 

These  data  represent  the  average  electron  density  along  the  zaxis. 

The  computation  of  both  spatially  and  temporally  resolved  electron  densities 
is  complicated  since  the  plasma  structure  has  motion  and  the  electromag* 
netic  fields  associated  with  it  are  time- varying. 

Inspection  of  Figures  V-15  to  18  reveals  that  the  trend  is  to  increase 

electron  density  as  voltage  is  increased.  At  500  microns  pressure,  how- 

15  3 

ever,  the  density  is  fairly  constant  at  about  6.  5  x  10  part/cm  regardless 
oi  voltage.  Overall,  there  seems  to  be  a  general  increase  in  electron  density 
as  time  goes  on.  There  is  one  exception,  though,  at  300  microns  there  is  an 
anomaly  between  12  and  l4kvfor  the  density.  At  other  pressures,  the  data 
aecms  more  consistant. 

One  other  independent  check  of  number  density  was  made.  The  line 

shapes  and  underlying  continuum  plots  made  for  calculating  temperature 

in  the  next  section  were  used  *o  measure  the  half  widths  of  H  and  H.  . 

Y  6 

See  Figures  V-19  and  20.  These  scans  were  made  in  a  manner  similar  to 

the  scan  shown  in  Figure  V-13  except  that  the  H  and  H  scans  are  to  a 
O  Y  6 

scale  of  3.  35 A/ inch  on  the  chart.  The  half  width  of  the  spectral  line  is 

measured  at  1/2  the  peak  line  intensity,  taking  care  to  correct  for  the 

underlying  continuum.  Shown  on  the  graphs  are  the  "0"  level,  vfoere  mono- 

chromater  shutter  is  closed,  and  the  level  of  continuum  under  the  line. 

This  last  level  is  a  computed  average  taken  from  the  far  wings  on  both 

sides  of  the  line.  The  calculated  half  width  for  H  and  H.  are  14.  2  A  and 

0  V  6 

12  A 


respectively. 
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The  electron  density  is  calculated  using  the  following  formula 
1  o 

from  Griem. 

N  =  C  (N  ,  T  )  U  %  5.  6) 

e  e  s 

where  AX  =  half  width  ir.  angstroms 

s  3 

N  =  electron  density  in  part/cm 

C  (N  ,  T)  a  coefficient  that  is  only  a  weak  function 
0 

of  density. 

14  14 

For  H  ,  CIN  ,T)*  3.  61  x  10  and  for  we  have  C(N  ,  T)  *  1  11x10 

y  e  o  e 


Therefore 


N  =  3.61x10  (14 

e 


14  (14.2)% 


N  =  1.  9x10 
e 


16  par V  3 

/cm 


cm  from  H 


N  «  1.11x10 
e 


U  (12)% 


N  =  4.7x10 
e 


15 -"/cm3 


cm  from  Hc 


^he  scans  from  which  these  densities  were  calculated  were  done  at  a 
sampling  time  of  2.0  microseconds  after  the  current  pulse  began.  The 
modulator  voltage  was  10  kv  and  the  neutral  pressure  was  500  microns. 

The  electron  density  measured  from  H  observations  for  these  conditions 

15  3  8 

was  6.7x10  part/ cm  . 


The  densities  computed  from  half  widths  of  and  are  consistant 
with  those  obtained  from  profile  analysis  of  Ha  jince  the  theory  fer  Stark 
broadening  of  and  H.  is  only  accurate  within  30%  while  that  for  is 
known  to  be  on  the  order  of  5%  .  The  merging  of  lines  at  essentially  the 

same  time  yielded  a  Inglis  -  Teller  calculation  for  electron  density  of 

15  3 

5.  7x  10  part/cm  . 
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Since  there  is  agreement  using  independent  checks,  the  measure¬ 
ments  of  electron  density  done  in  deta:'.  from  observations  of  Hg  should 
be  reliable.  One  last  check  was  mad  ,  however,  this  was  a  check  to  see 
if  the  plasma  was  optically  thin.  If  this  is  not  the  case,  then  self  absorption 
i3  taking  place  and  the  line  profiles  are  altered.  A  correction  must  be  then 
applied  in  order  to  calculate  the  correct  electron  density. 

Instead  of  the  usual  technique  of  placing  a  mirror  behind  the  plasma 

and  then  recording  doable  the  line  intensity  at  the  photomultiplier  output 

for  optical  thinness,  the  following  procedure  was  employe  A  strobe  lamp 

was  placed  behind  the  plasma  bottle  in  such  a  way  that  the  itrobe  lamp  light 

output  massed  through  the  plasma  and  into  the  monochromater  setup  .  The 

intensi. /  of  the  strobe  lamp  was  approximately  the  same  as  the  Hfi  output 

of  the  plasma.  The  strobe  lamp  output  for  a  1.  6 A  band  centered  at.  was 

recorded  separately  and  the  HR  output  of  the  plasma  was  recorded  separately 

o  c 

with  the  identical  1.  6  A  band.  Then,  both  sources  were  recorded  simultan¬ 
eously.  This  is  shown  as  Figure  V-21  ,  the  top  photograph  is  Ha  output,  the 

P 

middle  is  the  plasma  and  strobe  lamp  both  emitting  light,  and  the  bottom 
photograph  is  the  strobe  lamp  output  by  itself.  The  signals  shown  in  the  top 
and  bottom  photographs  when  added  together  do  result  in  the  signal  shown 
in  the  middle  one.  Therefore  the  plasma  is  optically  thin  for  Ha  . 

P 

C.  Measurement  of  Temperature 

A  time  resolved  spectrogram  of  the  Balmer  spectral  lines  and 
and  their  surrounding  continuum  was  made  using  a  motor  driven  monochro¬ 
mater  and  sampling  oscilloscope.  The  oscilloscope  output  was  fed  to  a  strip 
chart  recorder.  This  experimental  technique  has  been  described  in  greater 
detail  in  Chapter  IV-G. 

The  spectragrains  of  H  and  H.  are  given  in  Figures  V-19  and  V-20 

V  6  O 

respectively.  The  wavelength  scale  is  3.  35  A/inch,  and  the  spectrum  was 

taken  at  two  microseconds  after  the  beginning  of  the  current  pulse.  The 

o 

original  spectragrams  consisted  of  more  than  100  A  of  continuum  on  either 


» 
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side  of  the  spectral  line,  making  the  plot  several  feet  long.  In  order  to 
present  the  data  in  this  thesis,  the  line  wings  and  outlying  continuum  were 
sliced  off.  However,  the  figures  do  shown  the  average  value  of  continuum 
under  the  spectral  line  as  calculated  from  the  continuum  magnitude  100 /P 
away  from  the  line  core.  Also  shown  is  the  "zero"  level ;  the  level  on  the 
graph  which  corresponds  to  no  light  passing  through  the  monochromater. 

In  Griem's  book,  "Plasma  Spe-troscopy",  the  relative  line  to  contin¬ 
uum  ratio  is  given  on  page  279  as  the  following  equation  : 


3/7  7  7  f  m 

lL  3  "  rr  (137aQj  fgexp  <EH  -  EJ  /kT  , 

*C  ZXA\  g.  [(gff/2)  (kT/EH)  exp  (EH/n,2kT) 


+  ][(gf  J n*)  exp  (EH/n2kT)] 


where 


l  s  absorption  oscillator  strength  of  line 
g  =  statistical  weight  of  lower  state 
E^  ~  energy  of  lower  state 
Epj  *  ionization  potential  of  Hydrogen 

g..  =  statistical  weight  of  ion  ground  state 

g ^  ~  free-free  Gaunt  factor 
g^  =  free-bound  Gaunt  factor 


Also  n  is  summed  from  n  =  3  to  at  least  n  =  6  for  the  Baln.er  series.  Then 
t  20 

n‘=7  is  used,  according  to  Griem  .  The  magnitude  of  IT  /  I  versus 

JL«  C 

temperature  for  and  were  given  as  graphs  in  Figure  IV- 10  . 


The  spectrunis  shown  in  Figures  V-19  and  V-20  are  linear  plots. 
Therefore  I.  /I_  is  found  by  comparing  areas  under  the  curves.  This  is 
accomplished  with  a  planimeter.  First,  the  area  under  the  curve  extending 
w  50/C  from  the  spectral  line  core  was  measured.  This  represents  I  +I_ 
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/) 

Next,  areas  that  were  100A  wide  and  tangential  to  the  first  measured 
area  were  computed.  These  last  areas  were  averaged  and  labeled  as  I  . 
The  sketch  below  shows  the  situat  on. 


Application  of  the  above  procedure  for  the  spectragrams  of  and 

yield  I  /I  ratios  of  1.  38  and  .  363  respectively.  Using  Figure  IV- 10  we 
L  L» 

find  a  ratio  of  1.  38  for  corresponds  to  a  7.  5  ev  electron  temperature 
and  .  363  for  the  spectrum  yields  an  8.  5ev  electron  temperature.  The 
agreement  between  the  two  calculations  is  approximately  12%  . 


Now  we  have  to  check  to  see  if  these  temperatures  are  meaningful. 

Firs*  of  all,  are  the  measurements  in  error  due  to  a  failure  to  satisfy  the 

assumption  of  local  thermodynamic  equilibrium,  (LTE)  .  Well  according 
21 

to  Griem  ,  if  the  plasma  were  time-independent  end  homogeneous,  an 
electron  density  of 

18  1  /  kT  \  ^ 

Ne-7xl°  T 7L  ("IT!  )  is  required 
n  /2  H 


for  a  Hydrogen  plasma  so  that  LTE  is  satisfied  to  within  10%  .  For  a 

13 

temperature  of  8ev.  ,  then  the  electron  density  must  be  more  than  10  part/cm 
to  insure  LTE  within  10%.  Our  densities  are  on  the  order  of  5  x  lO^part/cm  . 
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therefore  the  validity  of  the  LTE  assumption  introduces  negligible  error 
for  a  plasma  structure  that  is  time- independent  and  homogeneous. 


Our  plasma  may  be  considered  to  be  transient  and  inhomogeneous. 

In  that  case  more  stringent  requirements  are  set  for  the  LTE  assumption. 
Treating  the  transient  case  first,  what  we  require  is  that  the  equilibrium 
time  for  satisfaction  of  LTE  be  about  10%  of  the  plasma  life  time.  In  our 
case  the  plasma  has  a  lifetime  on  the  order  of  10  microseconds,  therefore 
we  expect  that  the  LTE  equilibrium  time  must  be  less  than  0.  8  microseconds  . 


22 

The  following  equation  given  by  Griem  gives  the  equilibrium  time  for 
partial  LTE,  that  is,  for  all  principal  quantum  numbers  above  some  n  there 
can  be  total  LTE. 

7 

kT 

)  exP  l  „ 

n  kT 


T  R3 

n 


7  Hi 

4.5x10  /  kT 


2E 


4 

n  N 


(tt) 


sec 


Now  for  LTE  in  a  homogeneous  and  stationary  plasma  at  8ev  only  an 

N  =  10 13  was  required.  Substituting  these  values  into  the  equation  above 

6  -9 

for  all  n>5,  ,  etc.  ,  we  find  that  5.  7x10  sec.  The  conclusion 

is  that  deviations  from  LTE  between  excited  states  is  not  serious. 

Our  plasma  is  also  inhomogeneous.  In  addition  to  the  previous 

requirements  for  LTE,  we  now  must  show  that  the  gradient  of  electron 

23 

temperature  is  small  over  distances  that  a  particle  can  diffuse  in  times 
-8 

of  the  order  of  10  seconds  as  calculated  in  the  previous  section.  For  a 
temperature  of  8ev  the  electron  thermal  velocity,  V  is 


and  Ve  =  1.68x10  cm/sec.  Typical  diffusion  velocity  in  our  plasma  is 
lx  10^  cm/sec.  For  a  time  interval  of  1 0  ^  seconds  spatial  motion  due  to 
temperature  is  on  the  order  cf  1cm.  ,  while  that  due  to  diffusion  is  10  ^mm. 
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Since  the  magnitude  of  tne  temperature  spatial  gradient  must  be  inversly 
proportional  to  the  characteristic  path  length,  it  appears  that  the  spatial 
requirement  for  LTE  has  been  satisfied. 

All  in  all,  the  experimental  average  electron  temperature  of  8ev.  and 

15 

the  determination  of  electron  number  density  «  6x10  part/cm.  from  the 
previous  section  does  not  produce  a  paradox,  or  even  a  serious  refutation 
of  the  LTE  assumptions  that  are  prescribed  by  Griem  in  order  to  generate 
theoretical  working  equations  for  calculation  of  temperature. 


D.  Determination  of  Resistivity 

The  magnitude  of  electrical  resistivity  associated  with  the  plasma 

exerts  a  strong  influence  on  the  growth  rate  of  the  resistive  instability,  which 

has  been  observed  in  this  experiment.  The  resistivity  lias  been  determined  in 

two  ways.  First,  the  voltage-current  measurements  were  used.  The  initial 

voltage  across  the  modulator  capacitors  represents  an  electric  energy  source 
1  2 

containing  —  CV  Joules,  where  C  =  modulator  capacitance  and  V  =  initial 
voltage  across  the  capacitor  .  All  this  energy  is  assumed  to  go  into  ohmic 
heating  of  the  plasma  and  series  resistor,  R  ext.  The  energy  dissipated  in 
the  plasma  is  represented  by  an  equivalent  resistor,  having  the  value  <R>  . 
Knowing  the  approximate  geometry  of  the  plasma,  then  the  resistivity  can 
be  calculated,  given  <R>  . 


<n> 


A  <  R  > 
L 


5.  7) 


<n>  =  equivalent  resistivity. 

<R>  =  equivalent  resistance  . 

L  =  mean  circumference  of  plasma. 

A  =  cross-sectional  area  of  plasma. 

The  second  method  of  determining  resistivity  stems  from  the  temperature 

24 

measurements  made  in  the  previous  section.  In  this  case,  Spitzer  has 
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-% 

calculated  an  expression  where  r)  c  T  (T  is  electron  temperature)  and 
the  resistivity  is  only  a  weak  function  of  electron  density.  His  equation 
is  given  below. 


n  = 


5.21x10  Ln  A 


3/2 


where  A  =  1.55x10 


10  _T 
N 


3/2 


1/2 


5.8) 


q  =  resistivity  in  ohm-meters 
T  =  temperature  in  electron- voltB 
N^=  electron  density  in  particles/cm' 


When  the  experiment  was  being  conducted,  the  core  voltage  and 
primary  current  signals  were  always  displayed  on  a  Tektronix  555  dual-beam 
oscilloscope  .  These  displays  were  recorded  on  Polaroid  3000  ASA  film  for 
every  modulator  voltage  setting  and  also  for  the  various  primary  winding 
configurations.  Typical  voltage  and  current  data  are  shown  in  Figure  V-22, 
This  figure  shows  data  at  modulator  settings  of  8,  10,  12,  and  14  kilovolts  . 

The  neutral  gas  pressure  of  the  Hydrogen  was  500  microns  and  the  inner 
primary  winding  had  a  uniform  structure  so  that  the  current  carrying  plasma 
sheet  in  the  toroidal  bottle  broke  up  into  two  filaments. 

As  mentioned  previously,  to  calculate  resistivity  from  data  represented 
by  Figure  V-22  we  equate  electric  energy  of  the  modulator  capacitor  to  ohmic- 
heating.  This  relation  is  given  in  the  equation  below. 


j  CV2  =  (<R>  +  Rext)  j  I2(t)dt  5.9) 

C  -  modulator  capacitor 
V  -  peak  modulator  charging  voltage 
<R>  =  average  plasma  resistance 
Rext  -  external  modulator  charging  •’s’stor 
I(t)  =  plasma  current  as  function  of  tir.  ,  t. 
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The  value  of  C  is  .  0334  uf  and  of  R  ext  is  11.5  ohms .  However  when  R  ext  is 

reflected  from  the  primary  side  of  the  transformer  to  the  secondary  side, 

11.5 

the  plasma  region,  then  we  have  an  effective  R  ext  of  — or  .  068  ohns, 

n^ 

where  n  represents  a  transformer  turns  ratio  of  n=  13.  The  term  I(t)  is  the 
efective  plasma  current,  which  is  n  times  the  primary  current.  It  was 
observed  that  I(t)  is  very  nearly  equal  to  Iq  sinuut  for  0<oit<TT  according  to 
the  experimental  data.  Substituting  the  approximate  form  of  I(t)  into  equation 
5.  9)  and  integrating,  we  find  the  following  expression  for  <R>  . 

2 

<R>  =  ^  C  (  ”)  -  R  ext  5.  10) 

o 

Application  of  equation  5.  7)  results  in  an  expression  for  the  resistivity. 


<T|  > 


R  ext 


] 


or 


^  A  r  C  /  V  \  n  .  I 

<n>  *  117  J  '  Rext  j 


5.  li) 


where 


a  2 

A  rs  10  m 
L  3  x  10  *  m 

T  =  length  of  current  pulse  in  seconds  of  time  . 


The  graph  shown  in  Figure  V-23  illustrates  the  behavior  of  V  versus 
I  for  different  modulator  voltages  and  primary  winding  configurations. 
Examination  of  the  graph  indicated  that  the  ratio,  V/IQ,  decreases  for 
increasing  V.  According  to  equation  5.  11)  ,  we  come  to  the  conclusion  that 
<n>  ,  the  resistivity,  decreases  as  V  is  made  larger.  Offsetting  this  relat¬ 
ionship,  however,  is  the  dependence  of  T  upon  the  modulator  photographs. 
Inspection  of  Figure  V-22  reveale  that  T  decreases  for  increasing  V. 


In  any  case,  <T)>  was  calculated  from  equation  5.  11)  for  various 
V/I  ratios  and  their  corresponding  T  values.  The  results  are  shown  in 
Figure  V-24  where  V  is  the  ordinate  and  <tj>  is  the  abscissa.  An  independ¬ 
ent  check  for  these  values  of  <n  >  can  be  made  by  using  the  measured  values 
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of  temperature  and  electron  density  giv^"  in  sections  C  and  B  of  this 

chapter  respectively.  These  data  can  be  used  in  Spitzer's  formula, 

equation  5.  8) ,  in  order  to  calculate  <T)  >  .  With  a  modulator  voltage  of 

12  kv  and  neutral  gas  pressure  of  500  microns,  the  electron  density  was 
15  3 

found  to  be  6.  7  x  10  part/cm  and  the  temperature  due  to  line -continuum 
ratios  for  H  (  and  Ht  was  calculated  as  7.  5  and  8.  5  electron  volts  respect¬ 
ively. 


The  resistivity  according  to  the  above  data  was  found  to  be  the 
following. 


<T1> 


=  2.2x10 


ohm-meters 


H 

Y 


<n> 


1.8x10 


ohn-meters 


These  values  are  plotted  on  Figure  V-24  for  comparison  purposes. 

The  difference  between  the  average  resistivities  calculated  by  the  two 
techniques  is  less  than  20%  .  Inclusion  of  inductive  effects  for  equation  5.  9) 
would  decrease  the  calculated  resistivity  in  the  direction  of  that  determined 
from  temperature -density  measurements.  The  additional  complexity  for 
a  little  better  agreement  did  not  appear  justified  at  this  time  however  . 


E.  Observation  of  Rayleigh- Taylor  Instabilities 

In  this  section,  photographs  are  presented  that  indicate  the  onset  of  a 
Rayleigh- Taylor  instability.  The  neutral  gas  pressure  of  the  Hydrogen  gas 
in  the  toroid  was  50  microns  Hg.  The  photographs  were  taken  in  the  framing 
mode  with  the  camera  looking  into  ihe  side -wall  of  the  toroid  along  the  major 
axis. 

To  correlate  the  experimental  observations  with  the  theory  presented 
in  Chapter  II,  the  observed  wavelength  is  compared  with  the  predicted  wave¬ 
length.  The  predicted  wavelength  depends  on  the  plasma  acceleration  and 
the  viscosity  of  the  plasma.  Viscosity  depends  on  the  ion  density  and  temp- 
erture.  Since  the  system  is  in  LTE,  we  can  substitute  both  the  electron 
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density  and  temperature.  This  data  was  presented  in  sections  B  and  C  of 

this  chapter.  The  plasma  acceleration  has  been  measured  experimentally 

25 

and  compared  favorably  to  a  snowplow  model  by  K.  Stuart  .  Both  rectang¬ 
ular  and  cylindrical  geometries  were  used  as  theoretical  models  and  it  was 
found  that  the  difference  in  their  solutions  was  insignificant.  The  following 
is  a  derivation  of  the  snowplow  model  for  both  the  rectangular  and  cylindrical 
geometries. 

The  following  equation  is  for  a  snowplow  model  having  rectangular 
geometry.  A  sketch  of  the  geometry  is  shown  below 


The  model  assumes  that  all  plasma  is  swept  up  by  a  "magnetic  piston"  and  is 
confined  to  a  small  region  immediately  ahead  of  this  piston. 


where 


H  =  K 


5.  12) 


M  =  (m  +  m. }  n  x  =  p  x 
c  l  o  o 


Substituting,  we  have  the  following, 


5.  13) 
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vo 


The  following  is  the  snowplow  equation  for  the  cylindrical  case. 


^  _  _sL(  w  dr  N 
2  dt  \  M  dt  J 


5.  17) 


H  =  K 

M  =  p  F  (r2  -  R2) 
o 


Equation  5.  17)  then  becomes, 


d  /  _r 
dt  \ 


2-r2 


5.  18) 


The  corresponding  integral  equation  with  initial  condition  of  r  I  =R  is 

|t=o 

given  belcw  as 


Uf  /(,-£) 


5.  IV) 


By  letting  K(t)  =  ^sinuut  as  done  in  the  rectangular  model,  equation 
5.  19)  can  be  integrated  into  the  following  form. 


-  -  Lnf  j")  =  - -  (2n2t2+cos2|.t-  l) 

2  8,„20  R2  '  ' 

o 


5.  20) 


Since  this  is  a  transcendental  equation,  the  left  and  right  hand  sides  of 
equation  5.  20)  is  set  equal  to  some  function,  F. 


i.  e.  , 


(t)‘ 


-  Ln  (r/R) 


K 

~°  ° — r~  (  2uu2  t2  +  ccf  2'i't  -  1^ 

8oj  P  R 
o 
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where  K  =  10'"^n  ,  u.'=  10^  r/s  ,  p  =1.67x10  ^ 

o  o 

and  R  =  2.  5cm. 

These  two  parametric  equations  are  plotted  in  Figure  V-26.  Then  for  the 
same  value  of  F,  the  corresponding  T/r  and  t  can  be  found.  These  values 
are  plotted  in  Figure  V-25  and  form  the  locus  of  equation  5.20)  which  it, 
the  cylindrical  solution  to  the  snowplow  model. 

Figure  V-27  shows  a  streak  photograph  of  the  plasma  from  which 
the  experimental  parameter,  s  used  above  were  taken.  The  ordinate  of  the 
picture  is  radial  distance  and  the  abscissa  is  a  time  base.  The  inside  edge 
of  the  plasma  has  been  plotted  on  Figure  V-25  to  show  correspondence  with 
the  rectangular  and  cylindrical  models. 

The  conclusions  derived  from  Figure  V-25  are  two  fold.  First,  the 
difference  between  the  rectangular  and  cylindrical  solutions  are  insignificant. 
Therefore  use  of  the  rectangular  geometry  for  the  toroid  ie  justified.  This 
is  the  geometry  used  in  Chapter  II  for  the  derivation  of  the  instabilities. 
Second,  the  data  taken  from  the  streak  photograph  fit  the  model  very  well, 
and  so  the  experimental  plasma  must  be  undergoing  good  snow  plowing.  For 
early  times,  when  the  Rayleigh- Taylor  instability  is  developing,  the  radial 
trajectory  is  essentially  quadratic.  The  acceleration,  G,  in  this  region  is 
about  2  x  10^  ^s^  . 

The  sequence  of  framing  photographs  shown  in  Figure  V-28  show  the 
onset  of  Rayleigh- Taylor  instabilities  as  evidenced  by  the  scalloped  appear- 
ence  of  the  inner  edge  of  the  plasma  in  the  second  and  third  photos.  The 
exposure  time  for  each  photo  was  50  nano  seconds  and  the  time  interval 
between  photos  was  500  nano  seconds  .  Examination  of  these  photos  indicate 
a  wavelength  of  about  Xsslcm.  with  a  growth  rate  for  the  instability  less 
than  0.  5  microseconds. 

26  L  L 

v7e  can  use  Marshall's  expression  for  viscosity  which  ip  : 
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1.89xl015  Ti5^ 

V  s  .  1  1  ■  * . - 

n.  LnA 
i 

where  v  =  viscosity 

T.  =  T  =  temperature  in  electron  volts 
1  e  3 

n^  =  n^  =  density  in  particles /cm  . 

15  3 

Substituting  8  ev.  and  n.sa  5x10  part,  /cm  yields  Vw8.  Using 
equation  2.  40)  from  Chapter  II  fcr  the  predicted  minimum  wavelength  due 
to  viscosity  damped  Rayleigh- Taylor  instabilities  we  can  compare  with  the 
observed  wavelength. 

.  - 1/3  2/3  .10 

X min  =  2n  G  v  where  v«8,  G«s2xl0 

Solving  the  above  equation  we  get  that 

\  min  •=  .  93  cm  . 

This  is  in  good  agreement  with  a  value  of  1  cm.  found  in  the  experiment. 


F.  Observation  of  Resistive  Instabilities 

Initially,  experimental  investigation  of  the  toroidal  plasma  discharge 
was  directed  towards  the  radial  motion.  These  first  experiments  were 
performed  on  the  narrow  bottle  and  the  plasma  was  considered  to  be  uniform 
in  the  axial  direction  in  order  to  jimplify  the  construction  of  a  theoretical 
model.  The  validity  of  the  assumption  mentioned  above  was  checked  by  taking 
framing  pictures  of  the  plasma  :  i  the  r  -  z  plane. 

These  photographs  shown  in  Figure  V-29  show  a  time  sequence  of  0.5 
microsecond  intervals.  It  is  obvious  that  there  are  two  filaments  and  their 
separation  increases  as  time  progresses.  The  dark  vertical  bars  appearing 
in  the  photographs  are  the  individual  wires  making  up  the  outer  winding  on 
the  toroida..  bottle.  The  modulator  voltage  was  6kv  and  gas  pressure  was 
500  microns;  the  same  values  used  when  taking  magnetic  probe  measurements 
in  Chapter  V-A  .  Comparison  of  magnetic  flux  mappings  for  various  times 
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presented  in  Chapter  V-A  with  the  framing  photographs  of  plasma 
luminosity  prove  that  the  luminosity  coincides  with  the  location  of  the 
plasma  current  filaments. 

The  breakup  of  plasma  into  separate  rings  and  the  axial  behavior,  in 

general,  was  believed  to  be  a  type  of  resistive  instability.  In  order. to 

examine  this  phenomena  in  greater  detail  two  major  modifications  were 

made.  First,  the  theory  of  resistive  instabilities  assumes  a  uniform  sheet 

of  plasma  in  the  presence  of  both  a  perturbing  magnetic  field  and  either  a 

uniform  conductivity  or  a  spatially  varying  conductivity  across  the  thickness 

of  the  plasma  sheet.  It  is  desirable  to  have  an  experimental  plasma  that  is 

initially  uniform  in  keeping  with  the  theoretical  model.  However,  it  was 

noticed  that  the  photographs  i  Figure  V-29  did  not  really  appear  to  have  a 

uniform  configuration  at  the  beginning  of  the  t.me  sequence.  This  problem 

was  solved  by  the  introduction  of  r.f.  preionization  as  discussed  in  Chapter 

III.  The  preionization  of  the  toroidal  plasma  results  in  an  axially  uniform 

plasma.  Second,  only  two  ring/;  were  observed  when  conducting  experiments 

27 

on  the  narrow  toroidal  bottle  and  it  was  thought  desirable  to  excite  more 
than  two  rings  in  any  way  possible  so  that  the  data  might  be  more  compre¬ 
hensive.  To  this  end  a  second  bottle  was  constructed  that  was  twice  the  width 
of  the  narrow  bottle,  the  maximum  width  of  the  second  bottle  being  determined 
by  the  dimensions  of  the  transformer  core.  Also  it  was  found  that  inserting 
in  the  annulus  of  the  toroidal  bottle  a  primary  winding  with  spatial  periodicity, 
the  number  of  plasma  rings  that  developed  was  a  function  of  that  peri'  dicity. 
(The  details  of  primary  winding  construction  are  given  in  Chapter  ILL).  It 
was  found  that  if  the  primary  winding  is  divided  into  n  identical  coils  placed 
in  succession  along  the  axis,  then  the  plasma  would  split  into  n  +  1  separate 
rings  of  filaments. 

Experiments  conducted  on  thir  second  bottle  were  performed  at 
modulator  voltages  ranging  from  8kv  to  20kv  but  only  data  from  8kv  to  14kv 
was  compiled  since  higher  modulator  voltages  resulted  in  electrical  noise 
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so  large  that  the  data  was  deemed  unreliable.  Neutral  gas  pressure  was 
varied  between  1  torr  and  100  microns.  The  manifistation  of  resistive 
instabilities  as  filament  formation  was  evident  from  1  torr  to  about  300 
microns.  Below'  300  microns  radial  motion  becomes  dominent  and  at  about 
150-100  microns  evidence  of  Rayleigh- Taylor  instabilities  appear. 

The  data  that  follows  consists  of  photographs  depicting  development 
of  either  two,  three,  or  five  filaments  from  a  plasma  that  is  initally  axiallv 
uniform.  The  modulator  voltage  was  varied  between  6kv  and  14kv,  and  the 
neutral  gas  pressure  was  0.  5  torr  (500  microns)  . 

The  first  set  of  data  are  framing  pnotographs  of  a  plasma  developing 
into  two  filaments;  this  is  shown  as  Figure  V-30.  Photographs  showing  three 
and  five  filament  development  are  given  in  Figure  V-31  and  Figure  V-32 
respectively.  Modulator  voltage  was  12kvfor  the  plasma  photographs  just 
mentioned.  Photographs  were  taken  at  other  modulator  voltages  also. 

Of  prime  interest  is  the  quantitative  determination  of  growth  rate  for 
the  instability  that  has  been  photographed.  To  accomplish  this  aim  a  model 
has  been  devised  that  related  the  growth  rate  to  the  separation  between  the 
filaments  that  were  observed  on  the  framing  photographs. 

First  of  all,  the  theory  predicted  an  exponential  growth  rate,  e  vt 
with  uu>0,  for  the  instability.  Also  there  is  a  characteristic  wavelength 
associated  with  the  instability.  Using  the  above  information  the  plasma 
current  is  assumed  to  have  the  following  form. 

I(z,t)  =  I  -  i  (  e11*  -1)  cos  kz  5.21) 

o  o 

where  I  is  tha  plasma  current  in  the  absence  of  any  perturbation,  i^  is 
the  perturbation  current,  and  k  is  the  wavenumber  of  the  instability  . 


5.  25) 
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and 


1  f'F  _  o  xt  . 
—  ~  ="2x  —  e  cos  kz 


i  2 


l-(~)  f ext- 1  'jcos  kz 
\  I  /  »  / 


Now 

Therefore 


lim  ~  =  0  since  F(F,  t)  =  0  . 

-  dz 
z-? 

dz  >F/*t 

lim  —  =  ltm  - 

dt  *F  /  f'z 

z  «  5  z-*" 


or 


and  finally 


,  _ o  xt 

2(u  e  cos  kc 

d£  _  ■  o _ 

dt 


2k 


*o  (  eA*- 1 )  sink? 


•  -  IS 
?  ~  k 


o 

cot  k? 


We  can  integrate  equation  5.  28)  to  get  ?  =  ?(t)  . 


j  tankridr^  =  ^  f 
J  k  J  ,  -xt 


£o 


to  1  -e 


Equation  5.29)  integrates  to  the  expression  given  below. 


cos  kg 
cos  k£ 

o 


-LL’t 

-xt-  ln(l-e_Ujt) 

- 

-UL’t  -ln(l«e  °) 

o 

m  “ 

•  - 

If  t  is  the  minimum  time  for  5  to  be  realizable,  then  we  know  coskp  = 

°  xt  ° 

since  it  corresponds  to  R  =  I  -  i  (e  -1) 
r  oo 


In  cos  kP  =  -  ;ii(t- 1  )  -  In 
o 


1  -e 


-xt 


1  -  e"  iut° 


This  can  be  w  ritten  as 


cos  k- 


I  -R 

o 


xt 

where  e  -1 


I  -R 
_o _ 

i 

o 
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5.  27) 


5.28) 


5.  29) 


5.  30) 

1 


5.  31) 


7? 


or 


e-i 


-l 


cos 


I  -  R 
o 

1  (e  -1) 
.  o 


for  t>  t 
—  o 


5.  32) 


For  the  purposes  of  the  experiment,  we  examine  the  photographs 
and  determine  5^  and  5  at  times  t^  and  t^  respectively  .  A  more  useable 
form  of  equation  5.  32),  then,  is  a  difference  formula  that  yields  ^  ,  the 


growth  rate. 


uut,  ,  I  -R  . 

1,2  o  1 

e  =  — : -  - - 

l  cos  k^ 


+  1 


1,2 


or 


we 


I  -R 

out.  =  In  1 1  +  - 

1,2  l 


1 


cos 


1,2 


I  -R 

know  I  »R,  I  >i  and  therefore  -  ■ 11  >  1 
o  o  o  1 

o 


Also 


So 


1 


cos  k? 


>1  , 


1,2 

r  io-r 

**1,2  ”  lnl_r_ 


cos  k? 


1,2 
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cos  k§ 


and 


( It  « 


*2-*! 


In 


1 


coik?. 


5.  34 


The  equation  above  is  then  used  for  calculating  growth  rate  from  the  values  of 
5  and  t  measured  on  the  framing  photographs  .  The  time  interval  between 
framing  photographs  is  0.  5  microseconds,  sO  will  always  be  some 

integer  multiple  of  that  time.  The  values  of  ?  are  read  from  the  photo¬ 
graph  using  a  7  eye  loop  that  incorporates  a  scale  marked  off  in  0.  1  milli¬ 
meter  divisions.  The  width  of  the  bottle  is  12  centimeters  and  the  image  on 
the  photograph  is  28  millimeters,  therefore  all  readings  taken  with  the  eye 
loop  should  be  multiplied  by  '  l2%8  or  4.  3  times.  The  following  table,  Table  I, 
lists  the  values  of  ?  and  t  measured  from  photographs  taken  of  the  different 
filament  configurations  for  various  modulator  voltages. 
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Table  I 


Two  Filaments  (k  =  ,5cm 


kv 

ycm| 

AT(usec) 

u'x  10  ^(sec 

3 

1.25 

1.  75 

.  5 

4.8 

10 

1. 75 

2.25 

1.0 

3.8 

12 

1.5 

2.25 

1.0 

5.  3 

14 

1.0 

2.  0 

1.0 

4.8 

^(mm) 

Three  Filaments  ( k 

=  1. 045cm  * ) 

-5 

kv 

?2(mm) 

AT(usec) 

uix  10  (sec 

8 

7.53 

9.67 

.  5 

5.2 

10 

4.  1 

6.45 

.5 

2.  7 

12 

4.  3 

6.45 

.  5 

2.  8 

14 

5.  37 

7.  5 

.5 

3.6 

Five  Filaments 

(k  =  1. 6cm  1) 

kv 

^  (mm) 

?2(mm) 

AT  (usee) 

tux  10  (sec 

8 

5.  62 

7.  5 

.5 

1.  08 

lu 

3.  75 

5.  62 

.  5 

5.  7 

12 

5.  62 

7.  5 

1.0 

5.  4 

14 

3.  35 

5.  62 

.5 

6.  1 

Another  method  for  obtaining  the  growth  r"te  of  the  instability  is 
to  take  streak  photographs  of  the  plasma  instead  of  framing  photographs. 

Then,  a  microdensitometer  scan  of  the  streak  photograph  yields  the  vari¬ 
ation  of  light  intensity  generated  by  the  plasma  directly.  This  method  is 
discussed  in  Chapter  IV-E. 

Streak  photographs  of  two,  three,  and  five  filament  configurations 
were  obtained  for  various  modulator  voltages.  Examples  of  such  photographs 
are  given  in  Figures  V-33  ,  V-34,  and  V-35  for  two,  three,  and  five  filaments 
respectively.  The  modulator  voltage  in  these  last  three  figures  was  set  a; 

8kv  and  then  lOkvfor  each  of  the  figures. 

Shown  in  Figure  V-36  is  a  calibratici.  photograph  of  light  intensity 
obtained  with  a  photodiode  for  the  five  filament  case.  Simultaneously  a 
streak  photograph  was  taken.  The  microdensitometer  trace  for  this  last 
photograph  is  shown  in  Figure  V-37  and  finally  Figure  V-38  provides  the 
calibration  curve  that  is  used  to  translate  succeeding  microdensitometer 
scans  of  streak  photographs  to  equivalent  intensity  plots.  Separate  calibration 
curves  are  made  for  each  filament  configuration  and  in  particular  for  each 
experimental  run  to  insure  that  variations  in  exposure  time,  development 
time  of  film,  etc.  do  not  influence  the  data. 

A  typical  data  reduction  procedure  proceeds  as  follows.  First  the 
streak  photograph  is  marked  at  known  times  along  its  time  base.  Then 
scans  are  made  at  these  times  across  the  axis  of  the  image.  Typical  scans 
of  a  five  filament  streak  photograph  for  a  modulator  voltage  of  12kv  are  shown 
in  Figures  V-39  to  41  inclusive  and  represent  scans  one  microsecond  apart. 
Attention  is  now  fixed  on  one  of  the  filaments  and  its  peak  percentage  of 
transmission  and  minimum  value  of  percentage  transmission  between  two 
filaments  are  recorded  for  each  microdensitometer  scan.  These  values  are 
then  converted  to  relative  light  intensity  via  the  calibration  curve.  Next  the 
ratio  of  peak  intensity  to  minimum  intensity  is  formed.  A  ratio  equal  to  one, 


WBH 
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of  course,  indicates  growth  has  not  yet  begun.  Ratios  greater  than  on,:  are 
the  result  of  an  instability  growing.  Table  Li  below  illustrates  a  typical 
calculation. 


Table  II 


Five  Filaments  At  12kv 


%  Trans. 


Time(us) 

1 

2 

3 

peak 

3.6 

2.  0 

2.6 

min 

3.8 

2.2 

-  2 

Rel.  Arnpl,  Time(us) 

1 

2 

3 

peak 

2.37 

~1.4 

3.8 

min 

3.8 

-1.4  i 

i  2.75 

Time(us) 

i  ! 

2 

3 

ratio 

1.065 

~1.0 

1.  38 

Now  to  find  the  growth  rate,  which  is  supposedly  exponential,  the  ratio 

of  the  numbers  in  the  last  row  of  the  table  above  is  taken  for  the  times  of 

uut 

two  and  three  microseconds.  Then  1.38  =  e  where  t  is  one  microsecond. 

5  -1 

Solving,  the  growth  rate  must  be  id  =  3.  17x10  sec.  .  This  prescription 
is  followed  for  reducing  the  streak  photographs  for  every  filament  configur¬ 
ation  and  impressed  modulator  voltage.  A  summary  of  uch  calculations 
is  given  in  Table  III  as  a  listing  of  growth  rate  for  different  combinations. 
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Table  III 

Growth  Rate  Determined  From  Streak  Pictures 


Modulator  Voltage 

Two  F  il. 

Three  Fii. 

F  ive  F il, 

(kv) 

-5 

-5 

-5 

ju  x  1C 

ex  10 

d’X  10 

8 

1.  54 

1.  1 

1.  1 

9 

1.  55 

1.  07 

1.  7 

10 

1.  33 

0.  99 

1.  95 

11 

1.  35 

0.  95 

2.  5 

12 

1.  3 

0.  96 

3.  17 

13 

1. 75 

1.  1 

2.  51 

14 

1.6 

0.  98 

2.  87 

In  order  to  correlate  th?  observed  growth  rates  of  the  instability 
with  the  ther  reticai  behavior  of  resistive  instabilities  a  graph  is  plotted 
from  the  experimental  data.  The  graph  is  w  plot  of  p  versus  S  where  : 


and 


p  =  urR 

s  “  tr/th 


uoa 
tr~  <n> 


th= 


a(u  <P>) 

o 


Vz 


<B> 


Val  ues  of  iy  as  determined  from  framing  and  streak  photographs  are  given 
in  Tables  land  III  respectively.  The  resistive  diffusion  time,  t  ,  is  cal- 
culated  for  a  plasma  wioth  of  a  ^2  centimeters  and  various  <n>  that  were 
determined  in  section  D  of  this  chapter  for  different  modulator  voltages  and 
curren.  configurations.  Refer  to  Figure  V-24.  The  parameter  <  p>  =  rr..<n>, 
where  m.  is  the  Hydrogen  ion  mass,  has  been  evaluated  in  Chapter  V- B  for 
<n>  at  a  time  of  two  microseconds.  Since  this  is  approximately  the  same 
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instant  of  time  that  ju  was  evaluated,  calculations  of  <p>  are  acceptable. 
Values  of  <m>  are  selected  from  the  curve  for  500  microns  pressure. 
Magnetic  flux  density,  <B>  ,  is  calculated  from  equation  5.  5)  of  section 
A  ir  this  chapter.  For  convenience,  it  is  repeated  below. 

_  <I>  Weber /m^ 

<B>  *  %~nr 

2 

where  <  I >  -  —  lpeak 

TT  r 

L  =  equivalent  axial  length  of  plasma 

Table  IV  gives  a  tabulation  of  all  the  parameters  mentioned  above 

together  with  the  calculation  of  t_  ,  t  ,  S,  and  p  for  a  plasma  with  two 

R  H 

filaments  being  generated.  Similar  information  is  given  by  Tables  V  and 
VI  for  the  three  and  five  filament  case  respectively. 


Three  Filaments 


Mod.  Peak  Normalized  Growth  Rate 

Voltage  Current  Framing  Streak 


r 


8<* 


We  no’  have  p  versus  S  for  all  the  filament  configurations  and  these 

data  are  plotted  in  the  graph  shown  in  Figure  V-42.  The  data  taken  from 

framing  photographs  are  plotted  as  curves  distinct  from  those  reduced  from 

Streak  photographs  since  the  method  of  analysis  was  different  for  the  two 

2/5 

types  of  photos.  For  comparison  purposes,  a  plot  of  p  ■  S  is  included. 

This  corresponds  to  the  case  of  a  =  1  .  Also  included  is  a  curve  calculated 
28 

by  J.  Wesson  using  a  digital  computer.  He  set  a  =  1  and  did  a  numerical 
analysis,  making  no  approximations  lor  solving  the  basic  set  of  equations 


derived  by  F.  K.  R.  The  asymptotic  solution  found  by  F.K.R.  was,  of  course, 
+  2/5  2/5 

p  ~  a—  S  for  the  rippling  and  tearing  modes. 
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A. 

B.  Field  -  136  gauss/box 
radius  3.  656  cm. 

A. 

B.  Field  -  136  gauss/box 
radius  4.  1 56  cm. 

A. 

B.  Field  -  136  gauss/box 
radius  4.  656  cm 

A. 

B.  Field  -  68  gauss /box 
radius  5.156  cm. 

A. 

B.  Field  -  68  gauss/box 
radius  5.656  cm. 

A. 

B.  Field  -  68  gauss/box 
radius  6.  1 56  cm. 


time  scale  -  2usec.  /box  A.  Primary  current  145 amps,  /box 


B  Magnetic  Probe  -  Center  of  Discharge 
Figure  V  -  1 
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A. 


R.  Field  -  68  gauss/box 
radius  6.  656  cm. 


A. 


B.  Field  -  34  gauss/box 
radius  7.  1  56  cm. 


A. 


B.  Field  -  34  gauss /box 
radius  7.  5  cm. 


time  scale  -  2  usee,  /box 
A.  Primary  current  145  amps,  /box 


Figure  V-  l  Continued 
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A. 

B.  Field  -  32  gauss/box 
radius  3.  656  cm. 

A. 

B.  Field  -  32  gauss/box 
radius  4.  156  cm. 

A. 

B.  Field  -  32  gauss /box 
radius  4.  65^  cm. 

A. 

B.  Field  -  32  gauss/box 
radius  5.156  cm. 


A. 

B.  Field  -  32  gauss /box 
radius  5.  656  cm. 

A. 


time  scale  -  2  v.sec.  /box 


B.  Field  -  32  gauss /box 
radius  6.  156  cm. 

A.  Primary  current  14  Samps,  /box 


B  Magnetic  Probe  -  Center  of  Di.ch.rge 


Figure  V  -  2 
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A. 


B.  Field  -  16  gauss/box 
radius  6.656cm. 


A. 


B.  Field  -  16  gauss/box 
radius  7.  1 56  cm. 


A. 

B.  Field  -  8  gauss/box 
radius  7.  5  cm. 


time  scale  -  2  usee,  /box 
A.  Primary  current  145  amps,  /box 


Figure  V-2  Continued 


t 


91 


A. 

B.  Field  -  32  gauss/box 
radius  3.  656  cm. 


A. 


B.  Field  -  16  gauss/box 
radius  4. 1 56  cm. 


A. 


B.  Field  *  8  gauss/box 
radius  4. 656  cm. 


A. 

B.  Field  -  8  gauss /box 
radius  5. 1 56  cm. 


A. 


3.  Field  -  8  gauss /box 
radius  5.656  cm. 


A. 

B.  Field  -  8  gauss/box 
radius  6. 1 56  cm. 


time  scale  -  2  usee,  /box  A.  Primary  current  145  amps,  /box 


Bz  Magnetic  Probe  -  Center  of  Discharge 


Figure  V-  3 


B.  Field  -  Z  gauss/box 
radius  'j.  656  cm. 


Magnetic  Flax  Plot  -  Time  One  Microsecond 
Figure  V-4 


Magnetic  Flux  Plot  -  Time  Two  Microseconds 
Figure  V-5 


t 


F  'gure  V-6 


Magnetic  Flux  Plot  -  Time  Four  Microse 
Figure  V-7 
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Time  In  Microseconds 
0.0  1.5  3.0 

0.  5  2.  0  3.5 

1.  0  2.  5  4.0 


Framing  Photograph  of  Radial  Motion 


Figure  V  -  8 
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a.  SHEET  MODEL  b.  DUAL-RING 

MODEL 


Theoretical  Models  of  Current  Distribution 


Figure  V-9 
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current  -  4KA/box 
voltage  -  lKV/box 
time  -  1  usee,  /box 


narrow  slit 
ampl.  0. 5V/box 
time  Z  usee. /box 


wide  slit 
ampl.  5V/box 
time  Z  usee,  /box 


Photomultiplier  Tube  Output  Voltages 
Figure  V-  1 4 
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ri  0001 


MODULATOR  VOLTAGE 


g,_OI  X*N=  A1ISN3G  N0dl0313 


I4kv 


I 


200  300  400  500  600  700  800  900  1000 

PRESSURE  IN  MICRONS  OF  MERCURY 


i  10 


both  plasma  and  strobe 
light  radiating 

strobe  light  output  - 
1.  6  A  band  centered  at 


Measurement  of  Self  Absorption 


Figure  V  -  21 


Ill 


modulator  voltage 
8  KV 


10  KV 


Top  Trace  -  core  voltage  -  200V/cm 
Bottom  Trace  -  plasma  current  -  1250A/cm 
Time  scale  -  2  microseconds/cm 


Voltage  -  Current  Relationship 


F igure  V  -  22 
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Top  Trace  -  core  voltage  -  200V/cm. 
Bottom  Trace  -  plasma  current  -  125GA/cm. 
Time  scale  -  2  microseconds/cm. 


Figure  V  -  22  Continued 


2  FILAMENTS 

3  FILAMENTS 
5  FILAMENTS 
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KILO  -  AMPERES 


SNOWPLOW  MODEL 
FOR 

PLASMA  BOUNDARY 


-  30NV1SI0  1VIQVB 


Figure  V-25 
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Streak  Photograph  of  Radial  Motion 


Figure  V  -  27 


RADIAL  DISTANCE  -  cm 
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0.1  sec 


0. 6  n  sec 


l.i  f!  sec 


Framing  Photographs  of 
Rayleigh- Taylor  Instability 


F igure  V  -  28 
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Sequence  in 
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Narrow  Bottle  with  Outside  Winding 
(12KV  with  500u  pressure) 
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Figure  V  -  29 


outer  radius 


inner  radius 
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1.  5 

3.  0 

0.  5 

2.  0 

3.  3 

1.0 

2.  5 

4.  0 

Time  Sequence  in  Microseconds 


Two  Filament  Configuration 


Figure  V  -  30 
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outer  radius 


inner  radius 


0.  0 
0.  5 
1.0 


1.  * 
2.0 
2.  5 


3.0 

3.  5 

4.  0 


Time  Sequence  in  Microseconds 


Three  Filament  Configuration 


Figure  V  -  31 
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outer  radius 

inner  radius 


0.  0 

1.  5 

3.  0 

0.5 

2.  0 

3.  5 

i.  0 

2.  5 

4.  0 

Time  Sequence  in  Microseconds 

Five  Filament  Configuration 
F  igure  V  -  32 


cm 


Two  Filament  Streak  Photographs 
Figure  V  -  33 


Top  Trace : 
photodiode  output 


Bottom  Trace  : 

primary  current  150A/box 

Time  acale  2usec,/box 


Light  Intensity  via  Photodiode  Output 


Figure  V-  36 
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CALIBRATION  CURVE  FOR 
FILAMENTS 


Fl  ;’ure 


MiCRODENSITOMETER  SCAN  -  I  fjL  SEC 


— nzrnrzr:  nr.::.'  12 
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Figure  v-  k? 
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CHAPTER  VI 


CONCLUSIONS 


In  this  experiment  on  a  Hydrogen  toroidal  plasma  two  types  of 
instabilities  were  investigated. 

1.  A  Rayleigh- Taylor  instability  was  found  at  neutral  gas 

pressures  of  about  50  microns  of  Hg.  At  this  pressure  there 

1 W  2 

was  significant  radial  acceleration  of  the  plasma,  G«2xl0  s  , 
and  the  instabilities  were  observed  with  an  S.  T.  L.  image 
convertor  camera  looking  in  the  r  -  ♦  plane  . 

2.  Resistive  instabilities  were  observed  in  the  r  -  z  plane  as 
evidenced  by  an  axially  uniform  plasma  breaking  up  into  a  set 
of  filaments.  This  type  of  instability  was  found  for  neutral 
pressures  above  250  microns  Hg  and  detailed  measurements 
were  undertaken  at  a  pressure  of  500  microns  Hg. 

In  order  to  make  quantative  measurements  of  the  instabilities,  a 
number  of  plasma  parameters  were  measured.  In  short,  magnetic  field, 
electron  density  and  temperature,  induced  voltage,  plasma  current,  resist¬ 
ivity,  viscosity,  and  the  instablilty  growth  rate  were  found. 


The  magnetic  field  was  measured  both  inside  and  outside  the  plasma 
bottle.  Comparison  was  made  to  theoretical  models  with  moderate  success. 
In  any  event  the  prediction  for  average  magnetic  field  linking  the  plasma 
agreed  with  measured  values  to  an  error  of  about  5%  .  Electron  density  was 
measured  using  the  Inglis- Teller  method,  line  half- width  relations  utilizing 

H  and  H,,  and  a  sampling  of  the  line  profile  for  Ha  .  All  the  techniques 

^  “  15#/  1 5#/ 

yielded  electron  density  values  ranging  from  4.  7  x  10  /cc  to  6.7x10  /cc 

if  the  H  measurement  is  neglected  (this  gave  N  ~  1.9xl0^Vcc)  . 

Y  ® 
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Resistivity  was  calculated  from  V-I  probe  measuiements  and  also  from 
temperature  measurements.  Undor  the  same  operating  conditions,  those 
two  techniques  yielded  values  of  resistivity  within  20%  . 

Photographs  taken  of  Rayleigh- Taylor  instabilities  resulted  in  a 
wavelength  of  1cm.  compared  to  a  prediction  of  0.93cm.  which  is  about  10% 
in  error.  This  value  represents  good  agreement,  when  one  consideres  that 
an  idealized  linear  model  was  used  to  predict  phenomena  in  an  inherently 
non-linear  media,  i.e.  plasma.  Coupled  to  the  study  of  Rayleigh- Taylor 
Instabilities  was  the  study  of  radial  motion.  This  resulted  in  the  verification 
of  a  rectangular  snowplow  model  for  the  radial  trajectory  of  the  plasma 
found  by  experiment 

Observations  of  growth  behavior  of  resistive  instabilities  were  two 
fold.  One  set  of  growth  rate  data  and  macroscopic  behavior  was  compiled 
from  analysis  of  framing  photographs.  A  more  complete  set  of  data  was 
reduced  from  quantitative  streak  photography  using  a  microdensitometer. 
Resistive  instabilities  with  different  wavenumbers  were  excited  by  changing 
the  spatial-periodicity  of  the  primary  windings.  The  graph  shown  in  Figure 
V-  42  shows  all  data:  growth  rate  as  determined  by  framing  and  streak 
pictures  for  the  two,  three,  and  five  filament  configurations. 

The  curves  obtained  from  framing  photographs  are  grouped  much 
higher  than  those  obtained  by  streak  photograph  analysis.  The  relationship 
between  the  curves  for  different  filament  configurations  is  consistant 
between  the  two  photographic  techniques  however.  One  reason  that  the 
framing  data  indicated  larger  growth  rate  may  be  caused  by  an  insufficient 
gray  scale  for  Polaroid  film.  This  means  that  the  boundry  region  between 
plasma  and  vacuum  is  more  pronounced  and  makes  the  filaments  appear  to 
be  more  widely  separated  than  they  really  a:e.  Examination  of  such  a  film 
with  apparent  widely  separated  filaments  results  in  a  larger  calculated 
growth  rate.  On  the  other  hand,  data  reduction  of  streak  photographs 
employs  a  continuous  calibration  of  light  intensity  versus  film  density  so  the 
problem  mentioned  above  does  not  exist  for  this  case. 
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The  curves  reduced  from  streak  photographs  have  slopes  that  fall 
between  those  curves  due  to  F.K.R.  and  J.  Wesson.  In  particular,  the 
magnitude  of  the  experimental  curve  for  three  filaments  fits  between  the 
theoretical  curves.  The  theoretical  curves  were  for  a  normalized  wave- 
number  of  a  =  1  .  Interestingly,  from  Table  I  of  Chapter  V-F  we  see  that 
the  three  filament  case  corresponds  to  k  =  1.045cm.  ^  Since  the  plasma 

width  is  about  2cm.  ,  we  have  for  three  filaments  a  normalized  wavenumber 
of  order  a  ~  1  .  Values  of  k  for  other  filament  configurations  taken  from 
the  same  table  yield  normalized  wavenumbers  of  ct  ~  0.  5  and  a  ~  1.6  for 
two  and  five  filament  cases  respectively. 

According  to  the  analysis  of  Chapter  II  for  resistive  instabilities  the 
growth  rates  of  the  two  modes  considered  are  as  shown  below. 


Tearing  : 

P~a'2/5  S2/5  ;a<l 

Rippling: 

p  -  a2/5  S2/5  ;  a>)  ,  Vn  t  o 

A  particular  mode  can  therefore  be  identified  most  simply  by  inspection 

29 

of  the  inequality  that  a  must  satisfy  ,  Inspection  of  Figure  V-42  shows 
that  curves  for  two  and  five  filament  configurations  have  larger  growth 
rates  than  the  curve  for  three  filaments  ( a  ~  1  )  .  Since  a~  0.  5  <  1  for  the 
two  filament  case,  this  can  be  identified  as  mainly  a  tearing  mode.  Also, 
the  five  filament  case  has  ana-  1.  6>  1  and  is  therefore  associated  with  a 
rippling  mode  rather  than  a  tearing  mode  . 

Again,  the  linearized  theory  appears  adequate  to  describe  the  basically 
non-linear  behavior  of  a  resistive  instability.  The  magnitude  and  fu".  tijnal 
dependence  on  a  and  S  appear  to  be  satisfied. 
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u  abstract  visual  observations  i-nade  on  a  pulsed  toroidal  plasma  indicate  that  the 
current  carrying  plasma  sheet  breaks  up  into  filaments  parallel  to  the  direction  of 
the  current  and  perpendicular  to  the  direction  of  the  magnetic  field.  This  process 
involves  the  annihilation  and  fusion  of  initially  opposite  magnetic  flux  lines  and  leads 
to  a  final  configuration  which  corresponds  to  a  lower  energy  state.  Such  a  state, 
however,  would  not  be  accessible  in  the  absence  of  discipation  since  its  attainment 
requires  slippage  of  the  plasma  across  magnetic  flux  lines.  The  operating  mecha¬ 
nism  must  then  be  associated  with  a  resistive  instability  of  the  type  first  predicted 
by  Furth,  Kileen,  and  Rosenbluth  and  also  by  Coppi,  Greene  and  Johnson  in  Physics 
of  Fluids,  Volume  6,  1963.  Observations  of  resistive  instabilities  under  laboratory 
conditions  is  of  particular  interest  because  these  modes  are  held  responsible  for 
many  astrophysical  and  geophysical  phenomena  such  as  solar  flares  and  auroras  as 
well  as  anomalous  diffusion  rates  in  fusion  devices. 

The  experimental  setup  consists  of  a  plasma  filled  ring-shaped  container 
of  rectangular  cross  section,  where  a  discharge  is  induced  by  means  of  an  iron 
cored  transformer.  The  time  development  of  the  current  distribution  in  the  plasma 
is  determined  with  the  help  of  a  combination  of  magnetic  probes  and  an  S.T.L. 

Image  Converter  Camera.  Optical  spectroscopy  is  used  to  calculate  the  temporal 
variation  of  electron  number  density  and  temperature. 
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